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. Al volver a mi casa busqué el disco. No lo encontré. 

Hace años que sigo buscando." 
Jorge Luis Borges: El Disco. 
(El Libro de Arena, 1975.) 



i 



A los rios Vltava y Guatapurí: 
El encanto y lo sublime en el fluir eterno de dos aguas en un mismo Heráclito. 



ii 



A Manera de Gratitud 
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en sonrisas. En su contentura dice "¿les guardo almuerzo? quizá a que hora vendrán, no vayan a bebé 
mucho". Me cuentas sobre alguna nueva forma geométrica que ideaste para solucionar algún pro- 
blema de la Carpintería. En este momento se aguan los ojos de mi emoción que compartimos. Yo 
te cuento de mis desgastadas noches buscando solución a algún problema de mi trabajo. Y reímos 
y tomamos un trago de Old-Parr. Tu me insistes en que debemos agradecer siempre al pueblo y a la 
gente que nos ha ayudado a superarnos. Y encontramos amigos y se rebosan los vasos y se alegra 
la euforia. A este ritmo ya son las 9 de la noche y transcurrieron casi doce horas desde que salimos 
de la casa. Ya en ella, mi mamá nos sirve un chivo guisado y yuca cocida. En las conversaciones 
han quedado interminados muchos temas y tu opinión sobre quien será el próximo rey vallenato, casi 
nunca desaciertas. En este paso algún suceso me trae nuevamente al presente o a esta otra realidad o 
más bien a este otro sueño en el que estoy escribiendo los agradecimientos de mi trabajo y no puedo 
evitar la nostalgia y pienso que sigues conmigo, porque soy una extensión tuya. Mi nostalgia labró 
este reencuentro que ha sucedido muchas veces. 

Como debo seguir, expreso entonces mis gratitudes: 

A todo Borges y su poema "El Golem", que me mostraron Praga. A los múltiples libros que es las 
" Mil y Una Noches", A Cortázar y "Torito" que me han acompañado desde siempre. A Anthony 
Hopkins en "Hannibal", R.E.M y Louis Armstrong, por su continua cercanía a mi espíritu. En 
Praga: Al Vltava que trajo a mí a la Bella Teata dibujada en mágicas historias de Tahití y Praga. 
A cada rincón del antiguo Gheto, en el que busqué el Golem, uno de mis motivos de visita a la 
hermosa ciudad con su Kafka y sus fantasmas y Smetana. Al proyecto del Socialismo con Rostro 
Humano, que aunque disminuido por la opulencia Rusa del momento, me dejó ver la nostalgia y la 
fantasía Praguense. Al Grupo de Física Teórica de Charles University, que permitió injertar el mundo 
Europeo en mi provincia, en el sentido de J. Martí. Especialmente, a Mgr. Martin Scholtz, la pareja 
Mgr. Norman Gürlebeck and Julia, Mgr. Ivan Pshenichnyuk y Mgr. Michal Tarana, mis amigos en las 
muchas Cervezas, Vinos y Becherovka, y las pocas Física y Literatura que hemos podido compartir a 
pesar de mi intrincado english-costeñol. Al Prof. RNDr. Jifí Bicák, DrSc, dr.h.c. que me aceptó entre 
sus estudiantes y brindó nuevas luces. 

En Colombia: a cada una de las partículas que conforman los sistemas y los ríos. A mis amigos 
de la UIS, particularmente al GIRG, que han tolerado mi diaria mamadera e gallo. A la amistad, 
paciencia, humo, enseñanzas y "... tinto ..." desde el Profesor Guillermo González V. Como siem- 
pre, a la fuerza y desmedida amistad de Rafa, Zandra, Tito e Igno. A Laura R. y nuestros recuerdos 
de mandarinas y guayabas y múltiples lecturas. A la casa de los Nechos y sus malditas arepas con 
salchicha. También a los Nechos, Nando, Eduar, Jesu, David, Vivix, El Pelaito y Gerzzen. A Robbie 
Burns, Thomas Parr, Christopher, Goyo, Quiceno y Dixon, mis viejos hermanos. Al espíritu crítico 
del Profesor Oscar Neira. Ahora bien, la literatura me ha mostrado sus alas y me ha dado casi todo lo 
que tengo, incluyendo la Física. En ese proceso, me fueron mostradas las fraternidades de Hernando 
Motato, Luz Mary y Harold Paredes, a quienes expreso mis gratitudes. A Roberto y Armando. A la 
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literatura misma, con el maestro Eduardo Pastrana Rodríguez, que está en cada página que leo. 
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1.1. Naturaleza e Importancia del Problema 



a obtención e interpretación de soluciones de las ecuaciones de Einstein y de Einstein-Maxwell 



-/_ , . constituyen una tarea ardua e interesante. Durante los últimos años se ha obtenido, mediante 
diversos métodos, un gran número de soluciones de estas ecuaciones. Sin embargo, es valido afirmar 
que, un gran número de ellas no constituyen sino parte de un conjunto de "objetos de colección" 
dada su difícil interpretación o su limitada importancia física. Una de las características más funda- 
mentales de un sistema aislado en el universo es la simetría axial: soluciones de estas ecuaciones 
correspondientes a espaciotiempos estacionarios axialmente simétricos, por ejemplo, pueden ser úti- 
les para describir estrellas y galaxias en equilibrio termodinámico [59, 38]. Así pues, el estudio de 
soluciones, estáticas y soluciones estacionarias, de las ecuaciones de Einstein y de Einstein-Maxwell 
correspondientes a configuraciones discoidales de materia autogravitante axialmente simétrica es un 
problema realmente atractivo y que permanece abierto a la investigación. Presenta vasto interés en 
diversos aspectos relevantes tanto desde el punto de vista de la teoría de la relatividad en sí misma, 
como en el modelado y entendimiento de propiedades físicas de sistemas astrofísicos. 

De modo que, una gran suerte de esfuerzos ha sido dedicada al estudio teórico del campo gravita- 
cional generado por fuentes discoidales de materia que pueden modelar ciertos sistemas astrofísicos 
referenciados en la "literatura observacional". Los discos de materia sin presión, conocidos como dis- 
cos de polvo, por ejemplo, se estudian en astrofísica como modelos de cierto tipo de galaxias y discos 
de acreción. Del mismo modo, discos de polvo contrarrotante se estudian para modelar galaxias, por 
ejemplo 50 y 5 a [85]. El modelo de discos relativistas también presenta gran importancia en el es- 
tudio de superposición de discos de acreción alrededor de agujeros negros, y galaxias con agujeros 
negros supermasivos en su núcleo; objetos puramente relativistas [71, 48, 1]. De igual manera, dis- 
cos con presencia de campos electromagnéticos, especialmente campos magnéticos, juegan un papel 
importante en astrofísica en el contexto de estrellas de neutrones y binarias, enanas blancas y en la 
formación y evolución de galaxias. 

En el contexto de la teoría general de la relatividad, la construcción de soluciones autoconsisten- 
tes físicamente apropiadas para construcción de fuentes materiales y el campo gravitacional creado 
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por estas fuentes es un problema bastante complicado. Para encontrar tales soluciones, se tienen que 
imponer y solucionar el problema "externo", un problema con valor en la frontera impuesto sobre 
las ecuaciones de Einstein (ó de Einstein-Maxwell) cuya solución determina el campo externo, y el 
problema "interno" cuya solución determina la estructura y la dinámica de la fuente material en su 
propio campo gravitacional. Las soluciones físicamente razonables tienen que satisfacer algunas con- 
diciones adicionales tales como finitud y positividad de la masa, característica físicamente razonable 
de materia y extensión geométrica razonable. Existe un caso límite muy interesante. Con fuentes de 
discos infinitesimalmente delgados, el problema interno se reduce a la formulación de condiciones 
de frontera propias para la solución externa. Así, en este caso, soluciones autoconsistentes se pueden 
obtener solucionando el problema externo sujeto a algunas condiciones propias de frontera con el fin 
de obtener fuentes de materia con algún comportamiento físicamente razonable. 

Ahora bien, la obtención de soluciones de las ecuaciones de Einstein y de Einstein-Maxwell se 
puede llevar a cabo numéricamente, perturbativamente o de manera exacta. A pesar del progreso de 
la relatividad numérica y de que los computadores pueden manejar cierta variedad de asignaciones 
genéricas, las soluciones numéricas son difíciles de clasificar e interpretar. La solución analítica exac- 
ta, por su parte, representa un ideal que, no obstante, es solo asequible mediante simplificaciones 
sustanciales del problema. Sin embargo, la importancia de abordar problemas físicos, como los men- 
cionados en líneas anteriores, requiere cierta intuición que puede ser ofrecida por soluciones analíticas 
exactas particulares, las cuales a su vez pueden estimular el análisis de situaciones más generales, al 
tiempo que permiten discriminar más fácilmente entre una característica anómala y una característica 
física; así como también especificar dónde están las singularidades, cuál es su característica, si la so- 
lución es estable o no, o si es en algunos casos genérica. Asimismo, en el caso de la formulación de 
una conjetura respecto a una situación general inspirada en una solución analítica exacta particular, 
encontrar nuevas soluciones analíticas exactas puede jugar un papel importante, verificando, clarifi- 
cando o modificando la conjetura o, en algunos casos, las soluciones particulares pueden actuar como 
estados asintóticos de una clase general de modelos y llegar a ser aún más significativas. 

Los párrafos anteriores permiten ver la clara importancia que para el desarrollo de la teoría en 
sí y el conocimiento de las propiedades físicas de sistemas astrofísicos representa la obtención e 
interpretación de soluciones analíticas exactas de los sistemas de ecuaciones de Einstein y de Einstein- 
Maxwell en términos de modelos relativistas de discos delgados. 

1.2. Estado General del Tema 

Durante las últimas cuatro décadas la obtención y análisis físico de soluciones exactas estacio- 
narias axialmente simétricas de las ecuaciones de Einstein y de Einstein-Maxwell correspondientes 
a configuraciones discoidales de materia autogravitante ha suscitado una extensa variedad de traba- 
jos. Estos fueron estudiados inicialmente por Bonnor y Sackfield [8] y Morgan y Morgan [63, 62]. 
Posteriormente, diferentes autores han obtenido varios tipos de soluciones exactas correspondientes 
a discos delgados estáticos [93, 60, 13, 56, 50, 6, 7, 28, 22] y estacionarios [61, 5, 69, 29]. Se ha 
considerado, de igual manera, la superposición de un disco estático o un disco estacionario con un 
agujero negro [51, 52, 53, 76, 74, 77, 78, 42, 79, 80, 94, 96, 82, 81], discos de fluido perfecto con 
halo [90] , así como también se ha investigado la estabilidad de discos delgados que implican el uso 
de una perturbación de primer orden en el tensor de energía-momentum [88]. Fuentes discoidales 
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para espaciotiempos estacionarios axialmente simétricos con campos magnéticos han sido estudia- 
dos como fuentes para campos de Kerr-Newman [49], campos magnetostáticos axialmente simétricos 
[55], mientras que discos de fluidos perfectos cargados fueron estudiados por Vogt y Letelier [91], y 
discos de fluido perfecto cargado como fuentes de espaciotiempos tipo Taub-NUT por García-Reyes 
y González [19, 20]. 

En los casos anteriores, los discos se obtienen mediante un método tipo "problema inverso". Para 
discos sin campo electromagnético el método se puede describir de la siguiente manera: se toma una 
solución de las ecuaciones de Einstein, de modo que exista una discontinuidad en las derivadas del 
tensor métrico en la superficie del disco, posteriormente se obtiene el tensor energía-momentum a par- 
tir de las ecuaciones de Einstein. Finalmente, se estudian las propiedades físicas de la distribución de 
materia mediante un análisis del tensor energía-momentum así obtenido. El método se puede extender 
al caso en que los discos en consideración presenten campo electromagnético. En este caso, se toma 
una solución de las ecuaciones de Eistein-Maxwell en el vacío y, las condiciones exigidas sobre el 
tensor métrico en la superficie del disco son las mismas que para el caso anterior. Los discos también 
se pueden estudiar de manera "directa": se da el tensor energía-momentum de la fuente y luego se 
resuelven las ecuaciones de Einstein. Este último método ha sido utilizado por algunos autores para 
obtener discos de polvo [67, 44, 47, 46, 17, 45]. 

Ahora bien, mediante el método de problema inverso se puede no solo obtener, con relativa faci- 
lidad, discos delgados con tensores de energía-momentum que presentan suficiente generalidad, sino 
también usar esta construcción para lograr interpretaciones de soluciones de las ecuaciones de Eins- 
tein ( Einstein-Maxwell) conocidas, en el sentido que los discos pueden actuar como fuentes exactas 
para métricas de espacios dados por soluciones de vacío (electrovacío). De otro lado, cuando se usa 
el método de problema inverso en el caso estático, el tensor energía-momentum es diagonal y, su 
análisis es directo. Por ello, la mayoría de autores considera discos estáticos y, excepto para discos 
de polvo, todos los discos obtenidos poseen fuentes anisótropas con esfuerzo azimutal diferente del 
esfuerzo radial. De otro lado, cuando el espaciotiempo en consideración es estacionario, el tensor 
energía-momentum que se obtiene es no diagonal y el análisis de su contenido físico es menos trivial 
y, en general, la fuente obtenida es no solo anisótropa sino además con flujo de calor. Debido a este 
hecho, existen muy pocos trabajos acerca de discos estacionarios y, ellos están limitados, en general, a 
discos obtenidos con soluciones tales que este no presente flujo de calor [61,5, 69]. En realidad, solo 
en un trabajo se considera discos con flujo de calor, pero solo se hizo un análisis parcial del tensor 
energía-momentum [29]. 

También existe una interpretación alternativa de los discos obtenidos. Esta interpretación se logra 
mediante el modelo de contrarrotación (CRM), en el cual el tensor energía-momentum de la fuente se 
expresa como la superposición de dos fluidos perfectos contrarrotantes. Ahora, aunque esta interpre- 
tación se puede ver como meramente teórica, existe evidencia observacional de discos constituidos 
por corrientes de materia rotante y contrarrotante [73, 72, 4, 83]. En el caso estático, la interpretación 
del CRM no solo se puede usar para describir una fuente constituida por dos componentes isótropos, 
sino también para simular efectos rotacionales verdaderos. Un estudio detallado del CRM para discos 
delgados estáticos relativistas genéricos se presenta en [22], aquí los autores obtienen expresiones 
para todas las cantidades físicas que caracterizan el modelo. Hasta ahora, sólo se ha logrado cierto 
avance en lo concerniente al estudio del CRM para fuentes de campos electromagnéticos y gravita- 
cionales que actúan como discos sin presión radial [18] y discos con presión radial [20]. En ambos 
casos los discos son estáticos y utilizan para su interpretación física, el método de problema inverso. 
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De otra parte, el CRM para discos delgados estacionarios no se ha desarrollado completamente y 
solo una versión preliminar de éste se presenta en la referencia [29]. En el trabajo mencionado, los 
autores consideran que el CRM se puede aplicar solo en casos en los que no existe flujo de calor y 
el esfuerzo azimutal es positivo (presión). Estas condiciones, sin embargo, son muy restringentes y, 
en muchos casos, se tienen discos que concuerdan con ello solo en una región parcial. Entonces, el 
CRM será válido solo como una interpretación parcial de los discos correspondientes. En un trabajo 
en desarrollo [26] los autores demuestran que estas condiciones no son necesarias e interpretan el 
CRM para discos delgados estacionarios axialmente simétricos con flujo de calor y esfuerzo azimutal 
arbitrario. 

Ahora bien, resultan particularmente interesantes los trabajos correspondientes a discos de polvo 
en espaciotiempos conformestáticos [92, 43]. En este caso, la densidad de carga del disco es igual a 
su densidad de masa y así, las fuerzas eléctricas y gravitacionales se encuentran en perfecto balance. 
Esta clase de configuraciones de equilibrio han sido llamadas por algunos autores "discos eléctrica- 
mente contrapuestos" (ECD) y han sido estudiados con cierto detalle tanto por teorías clásicas como 
relativistas en [1 1, 9, 10, 32, 89]. Ahora, como el contenido material de la fuente es polvo, probable- 
mente el tensor energía-momentum correspondiente estará siempre en acuerdo con las condiciones de 
energía, un hecho que es particularmente importante para los modelos de discos relativistas delgados. 
En efecto, como se puede ver en algunos de los trabajos anteriormente mencionados, existen muchos 
modelos de discos relativistas delgados que no están de acuerdo con estas condiciones y también mu- 
chos que cumplen sólo parcialmente con ellas. Los discos delgados conformestáticos presentados en 
[92, 43] fueron obtenidos por medio del bien conocido método de "desplazamiento, corte y reflexión" 
con el fin de introducir una discontinuidad en la primera derivada de una solución suave. El resulta- 
do será entonces una solución con una singularidad tipo función delta de Dirac en la hipersuperficie 
z = 0, que puede ser inteipretada como un disco infinito delgado. De otro lado, las soluciones que 
pueden ser interpretadas como discos delgados de extensión finita se pueden obtener si se introduce 
un sistema coordenado propio. Un sistema coordenado que se adapte naturalmente a una fuente fini- 
ta y presenta las condiciones de discontinuidad requeridas son las coordenadas esferoidales oblatas. 
Algunos ejemplos de soluciones que representan discos delgados finitos expresados en términos de 
estas coordenadas se pueden encontrar en [8, 63, 93, 56]. 

Los datos observacionales que soportan la existencia de agujeros negros en los núcleos de algunas 
galaxias, incluyendo los que registran la fuerte dinámica asociada a la evidencia dinámica generada 
por la fuerte interacción del resto de la Vía Láctea con su centro, de modo que no hay duda de la rele- 
vancia del estudio de sistemas binarios compuestos por un disco delgado que rodea un agujero negro 
(ver [12, 3] si se quiere apreciar recientes revisiones sobre datos observacionales). Por consecuen- 
cia, durante los últimos años, un gran número de trabajos ha sido desarrollado con el fin de obtener 
cada vez mejor información acerca de los aspectos relacionados con la dinámica de estos sistemas 
(ver [75, 42] para una revisión de los trabajos más relevantes). Ahora, debido a la presencia de un 
agujero negro, los campos gravitacionales presentes son tan fuertes que el marco teórico propio para 
el estudio de estas configuraciones lo constituye la teoría general de la relatividad. Por consecuencia, 
un gran esfuerzo ha sido dedicado a la obtención de soluciones de las ecuaciones de Eistein y de 
Einstein-Maxwell correspondientes a fuentes discoidales delgadas con un agujero negro central. 

Las soluciones estacionarias y axialmente simétricas de las ecuaciones de Einstein y de Einstein- 
Maxwell constituyen la mejor elección cuando el fin es describir analíticamente los campos gravita- 
cionales debido a discos alrededor de agujeros negros. Asimismo, tales espaciotiempos son de clara 
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importancia astrofísica, dado que ellos describen la configuración de equilibrio en la región exterior 
de sistemas rotantes. Así que, a través de los años muchos ejemplos de soluciones correspondientes a 
agujeros negros o a fuentes discoidales delgadas han sido obtenidos mediante diferentes técnicas. Sin 
embargo, debido al carácter no lineal de las ecuaciones de Einstein, soluciones correspondientes a la 
superposición de agujeros negros y discos delgados no son fáciles de obtener. De modo que, hasta 
ahora, soluciones exactas estacionarias correspondientes a este tipo de configuraciones no han sido 
obtenidas. 

Si se consideran discos delgados que se extienden más allá del horizonte de eventos, la materia 
localizada cerca el agujero se moverá con velocidades superlumínicas, tal como fue demostrado por 
Lemos y Letelier [51, 52, 53]. Así, con el fin de evitar la aparición de materia tachyonica, los discos 
delgados deben de tener un borde interno con un radio mayor que el radio fotónico del agujero negro. 
Así que, el problema de valor en la frontera para la ecuación de Laplace es más complicado y por 
consecuencia muy pocas soluciones han sido obtenidas. Esta clase de soluciones fue estudiada pri- 
meramente por Lemos y Letelier [52] haciendo una transformación de Kelvin con el fin de invertir la 
familia de discos finitos de Morgan-Morgan [63]. Ahora, aunque la segunda función métrica de esta 
solución no puede obtenerse analíticamente, su principal propiedad fue analizada ampliamente en una 
serie de trabajos por Semerák, Zácek y Zellerin [82, 81, 80, 79, 74, 94, 77], mediante computación 
numérica. 

Aparte de los discos invertidos de Lemos y Letelier, sólo otras dos soluciones para discos estáticos 
con borde interno han sido obtenidas. La una con discos isócronos invertidos [44] y la otra con una 
"power-law density" [78]. Mientras que Zellerin y Semerák [96] obtuvieron una superposición esta- 
cionaria usando el método de Inverse-scattering de Belinskii-Zakharov; no obstante, el análisis de sus 
propiedades es complicado dado que las funciones métricas no pueden ser obtenidas analíticamente. 
Además, esta solución implica una superficie sin sentido físico definida entre el horizonte del agu- 
jero negro y el disco [75]. Finalmente, una clase general de soluciones estacionarias fue presentada 
por Klein [46], mediante el uso de técnicas de superficie de Riemann, por medio de la cual sistemas 
disco-agujero negro físicamente aceptable se puede hallar, en principio. 

Ahora, una característica común entre todas las soluciones anteriormente mencionadas es que sus 
funciones métricas no pueden ser completamente calculadas analíticamente. Entonces, el análisis de 
sus propiedades matemáticas y físicas es muy complicado. Además, casi todas las soluciones presen- 
tan singularidades en el borde del disco, quizá la única excepción la constituye el tipo de soluciones 
presentadas por Klein [46]. La presencia de estas singularidades en casi todas las soluciones obtenidas 
ha sido considerada por algunos autores estrechamente conectadas con la no física delgadez infinita 
de la fuente. Sin embargo, tal como se verá más adelante, si el correspondiente problema de valor en 
la frontera se soluciona apropiadamente, es posible obtener una familia completa de soluciones libres 
de singularidades. 

1.3. Propósito y Organización del Presente Trabajo de Investiga- 
ción 

El objetivo central del presente trabajo es la construcción de modelos relativistas de discos del- 
gados a partir de soluciones exactas de las ecuaciones de campo de Einstein y de Einstein-Maxwell. 
Para tal fin, en primera instancia, se considerarán soluciones axialmente simétricas de dichas ecua- 
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ciones en las cuales la primera derivada del tensor métrico y del potencial electromagnético presentan 
una discontinuidad finita a través de un disco delgado. Dada entonces una solución con las propie- 
dades anteriores, se obtendrán el tensor de energía-momentum superficial y la densidad de corriente 
superficial de la distribución de materia por medio del formalismo necesario para el tratamiento de 
las ecuaciones de Einstein y de Einstein-Maxwell en términos de distribuciones tensoriales. 

Las propiedades físicas de la distribución de materia se estudiarán entonces mediante un análisis 
detallado del tensor de energía-momentum superficial, considerando soluciones suficientemente ge- 
nerales de las ecuaciones de Einstein-Maxwell que describan una distribución superficial de materia 
no sólo anisótropa, sino que además presente flujo de calor. Igualmente, se presentará una interpre- 
tación alternativa de tales discos con base en el modelo de contrarrotación, en el cual el tensor de 
energía-momentum superficial y la densidad superficial de corriente se expresan en términos de la 
superposición de dos fluidos perfectos contrarrotantes. 

Así, el Capítulo 2 expone un formalismo suficientemente general a partir del cual se interpre- 
tan ciertas soluciones de las ecuaciones de Einstein-Maxwell en términos de modelos relativistas de 
discos delgados. Se presenta un análisis detallado del tensor energía-momentum superficial (SEMT) 
correspondiente discos delgado estacionarios axialmente simétricos. Se considera la clase mas gene- 
ral de discos delgados estacionarios con presión radial y flujo de calor diferentes de cero. Se analiza 
el contenido físico del SEMT y se obtienen las expresiones para la velocidad, densidad de energía, 
densidad de corriente eléctrica, esfuerzos principales y flujo de calor. También se presenta el Mode- 
lo de Cotrarrotación (CRM) con el fin de interpretar estos discos considerándose el SEMT como la 
superposición de dos fluidos perfectos contrarrotantes. Se obtiene una condición sobre la velocidad 
tangencial, necesaria para la descomposición del tensor energía-momentum, y se obtienen entonces 
expresiones para las densidades de energía y las presiones. Se considera la posibilidad de contrarrota- 
ción a lo largo de electrogeodésicas, así como la contrarrotación con velocidades iguales y opuestas. 
Mediante la escogencia dada para las velocidades de contrarrotación, se obtienen expresiones explí- 
citas para las densidades de energía y presiones de los fluidos contrarrotantes. 

En el Capítulo 3 se presenta una familia infinita de discos de polvo cargados eléctricamente, axial- 
mente simétircos de extensión finita. Con el fin de obtener estas soluciones de discos delgados finitos, 
se solucionan primero las ecuaciones de Einstein-Maxwell para espaciotiempos conformestáticos asu- 
miéndose que la función métrica y el potencial eléctrico mantienen una dependencia funcional y que 
la función métrica depende funcionalmente de otra función auxiliar, la cual se toma del conjunto de 
soluciones de la ecuación de Laplace. Como se puede ver en [33] este procedimiento es bastante útil 
dado que se pueden obtener "nuevas soluciones" a partir de ciertas soluciones "semilla" apropiadas y, 
el álgebra en el proceso es particularmente simple. Del conjunto soluciones de la ecuación de Laplace 
conocido, se toma aquí como función auxiliar las dadas por la familia infinita de discos de Kalnajs 
generalizados obtenidos por González y Reina [30], expresadas en términos de las coordenadas esfe- 
roidales oblatas y que corresponden a una bien comportada familia newtoniana de discos axialmente 
simétricos de radio finito. 

En el Capítulo 4 se presenta una familia infinita de soluciones exactas nuevas correspondientes a 
discos delgados con un borde interno central constituidos por polvo. Estas soluciones describen discos 
cuya densidad de energía es siempre positiva y bien comportada, de modo que su tensor energía- 
momentum está en completo acuerdo con todas las condiciones de energía. Ahora, aunque los discos 
son de extensión infinita, todos ellos tienen masa finita. Además, las soluciones son asintóticamente 
planas y regulares en todas partes, tal como se demostrará mediante el cálculo de todos los escalares 
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de curvatura. En dicho capítulo también se presenta la superposición completa de esta familia de 
soluciones con un agujero negro de Schwarschild, ampliando los detalles presentados por González 
and Gutiérrez-Piñeres en [31]. La solución obtenida representa un disco anular delgado infinito (un 
disco con un borde interno) alrededor de un agujero negro de Scchwarzschild. La masa del disco 
es finita y el tensor de energía-momentum cumple con todas las condiciones de energía. La solución 
puede interpretarse también como si describiera un disco delgado constituido por dos fluidos de polvo 
contrarrotantes que satisfacen también las condiciones de energía. 

En el Capítulo 5 se presenta, lo que constituye, hasta donde se sabe, la primera familia infinita de 
soluciones exactas completamente integradas explícitamente de las ecuaciones de Einstein-Maxwell 
correspondientes a la superposición de un disco de polvo contrarrotante con un agujero negro. Las so- 
luciones obtenidas representan la superposición cargada de un disco infinito delgado anular alrededor 
de un agujero negro de Schwarzschild. La masa total del disco es finita y el tensor energía-momentum 
está de acuerdo con todas las condiciones de energía. Además los estados asintóticos de la superpo- 
sición cargada se reducen a la superposición sin carga eléctrica descrita en el capítulo anterior. La 
solución puede interpretarse como la descripción de un disco constituido por fluidos provistos de car- 
ga eléctrica contrarrotando alrededor de un agujero negro y que, están también en completo acuerdo 
con las condiciones de energía. Como parte final se presentan las principales conclusiones del trabajo. 
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2.1. Introducción 

Con el fin de obtener soluciones de las ecuaciones de Einstein-Maxwell interpretables como dis- 
cos delgados estacionarios axialmente simétricos con campo electromagnético se puede tomar, 
como punto de partida, el elemento de linea de Weyl-Lewis-Papapetrou [84], en coordenadas cilin- 
dricas xf = (í, ip, r, z), 

ds 2 = - e m {dt + Wdcpf + e- 2 °P? 2 íV 2 + e 2A (dr 2 + dz 2 )], (2.1) 

donde O, A, "W y R son funciones de r y z únicamente. Por naturaleza cilindrica de las coordenadas 
[63] se entiende que r = 0 sobre el eje de simetría y, para z fijo, r crece monótonamente al infinito; 
mientras que, para r fijo, z crece monótonamente en el intervalo (-oo, oo). La coordenada ¡p varía en 
el intervalo usual [0, 2n). 

El sistema de ecuaciones de Einstein-Maxwell en el vacío, en unidades geometrizadas tales que 
c = &7tG = [i 0 = e 0 = 1, está dado por 

G" b = T" b , (2.2a) 
T llb = F m F b c - l -g ah F, t ,F"', (2.2b) 
F ab b = 0, (2.2c) 
F u b = d a A„ - d b A„, (2.2d) 



Capítulo 2. Modelos Relativistas de Discos Delgados con Campo Electromagnético 



donde todos los símbolos tienen el sentido usual. 

Para la métrica (2. 1), el sistema anterior puede ser escrito explícitamente como 

■K , + 7? = 0, (2.3a) 

e 44, 7T 2 [7?('H'A. r - 'W 1 f,r\r + 7?CH>A Z - W 2 ^-)- - CWA, r - WV,,W, - 27?$,.) 

- (<WA¿ - 'W i'/ )(7? - 27?<D, ; )] + (Ki//, r \ r + («(/'.-)„- - 2«(<A. r O,, + á Ú> ) = 0, 

(2.3b) 

7? [( r Wip. r ), r + CW^h] + (7?,, - 27?<í>, r )(A, r - <Wi//s) + (7?, z - 2«<D, ; )(A, ; - 

- 7?(A_, T + A „) = 0, (2.3c) 

e m (<W\ + IV 2 ) - e 44> [(A.,. - r W^, r ) 2 + (A¿ - WiftJ 2 ] + 2He m [(7?<D.,.),,. + CK'l> ) | 

- 7? 2 (<A 2 + iA 2 ) = 0, (2.3d) 

e 2 * [4(o,,nv, r + + mn^'W,), + KffrVjJ - 2(A^, + a^j 

+ 21V((/r + iA 2 ) = 0, (2.3e) 

e «»(^2 _ <w 2 ) _ 2e 2 * [(A,,. - r W\ii, r f - (A, - W>f,J 2 ] + 27?V 2Í '(^ - <A 2 ) 

+ 47? [(A,r7?,r - A V? ) - 7?(<D 2 - í> 2 ) + 7? J = 0, (2.3f) 

e 4 "^/^ - 2<r" > (A., - nViA,,.)(A., - -W<frJ + 27?[A, r 7?, z + A,-7?. ,. - 7? 

- 27?0,4>.-] + 27?V M ^ ir ^ = 0, (2.3g) 

donde A„ = (0-, A, 0, 0) es el 4-potencial electromagnético, cuya dependencia funcional es sólo de las 
variables ryz- 

Ahora bien: de (2.3a) se puede ver que la función 7?(r, z) satisface la ecuación de Laplace, de modo 
que se la puede pensar como la parte real de una función analítica T(£) = 7?(r, z) + iZ(r, z), donde 
f = r+iz. Por consecuencia, la función f(( ) define una transformación conforme de las coordenadas, 

r 7?(r,z), 

(2.4) 

z -> Z(r, z), 

de tal forma que la métrica (2. 1) puede ser escrita, en las nuevas coordenadas, como 

ds 2 = - e w {dt + !Rd<pf + e-^P? 2 V + e 2n (d-R 2 + dZ 2 )l (2.5) 
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donde el conjunto de funciones 





= Rer(f), 


(2.6a) 


Z(r,z) 


= Imr(f), 


(2.6b) 


1>0, z) 


= ¥(7?,Z), 


(2.6c) 




= íl(?i,Z), 


(2.6d) 


A(r,z) 


= n(7?,Z) + inF'(a 


(2.6e) 


ifrir, z) 


= «3, 


(2.6f) 


A(r,z) 


= A(K,Z), 


(2.6g) 



satisfacen el sistema 

+ 2K<F Z C?IA, Z - ^V,z)] + (^),« + WzXz " + íAz^z) = 0, (2.7a) 

■R- [ ^, K - W,«),»¡ - Wz)z + «(« _1 A, K ),« + ftift-'A^lz 

+ 2>F,«(A,« - 3\>fj^) + 2>I'. Z (A, Z - jTIi/^) = 0, (2.7b) 

é> 6 W K + ^I 2 Z ) - [(A« - ^, K ) 2 + (A z - ^(A, z ) 2 ] - K 2 ^ + tf z ) 

+ IKe 2 ^ [WVjdx + (Wz)z] = 0, (2.7c) 

e 2V [4(>P,^, K + f¿n x ) + K(K- ] ^, K ):r + W^zXz] + 2^(<A 2 R + 1?£> 

-2(A. K ^ + A, ZÍ A. Z ) = 0, (2.7d) 

- ^ 2 Z ) - 2e 4ip [(A« - ¿ñif/^) 2 + (A, z - ^(A, z ) 2 ] + 2ft 2 ((A 2 K - iA 2 z ) 

+ Afte™ [n, R - KC¥% - W 2 Z )] = 0, (2.7e) 

e m 3\^3\, z - 2e^(A :R - ^, K ){A X - 3\>fj x ) + le^ft [n, z - l»^] 

+ 27?V, K <A, Z = 0, (2.7f) 

que representa el conjunto de ecuaciones de Einstein-Maxwell en el vacío para un espaciotiempo 
estacionario axialmente simétrico en coordenadas canónicas de Weyl, cuyos métodos de solución se 
conocen ampliamente [84]. 

Ahora, con el fin de obtener una solución interpretable como el campo exterior de un disco del- 
gado con campo electromagnético, se supone que O, A, "W y R son funciones simétricas de z y que 
las primeras derivadas del tensor métrico son discontinuas en el plano z = 0. Soluciones con este 
tipo de comportamiento se pueden obtener por diferentes caminos mediante la elección apropiada de 
soluciones del sistema (2.7) y la función analítica T(l,)- 

Un procedimiento comúnmente utilizado es el denominado método de "desplazamiento, corte y 
reflexión". El método consiste en lo siguiente: se toman f(0 = f y cualquier solución arbitraria 
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del sistema (2.7), la cual en general no posee la simetría de reflexión requerida. Luego, con el fin de 
obtener la simetría de reflexión requerida, se hace la transformación z — > |z| + Zo, donde zo es una 
constante positiva. El resultado será entonces una solución con una singularidad tipo función delta 
en z = 0, la cual se puede interpretar como un disco delgado infinito con campo electromagnético. 
Algunos ejemplos se pueden ver en las referencias [18, 55, 91]. Mientras que soluciones interpretables 
como discos delgados infinitos pero sin campo electromagnético, obtenidas mediante este mismo 
procedimiento, se presentan en [5, 6, 7, 29]. 

Las soluciones logradas mediante el método de "desplazamiento, corte y reflexión" son, en cierto 
sentido, artificiales; dado que estas se obtienen por medio de una transformación geométrica. Una 
clase de soluciones "naturales" se puede obtener apropiadamente mediante la solución de un proble- 
ma de frontera. Para ello, se toma T(£) = ( y se introducen las coordenadas esferoidales oblatas, las 
cuales son adaptadas naturalmente a una fuente discoidal y cumplen con el requisito de la discontinui- 
dad. Posteriormente se escogen soluciones del sistema (2.7) que satisfagan las condiciones de frontera 
apropiadas en coordenadas esferoidales oblatas. Las soluciones obtenidas se pueden interpretar como 
discos delgados finitos sin presión radial y, en ausencia de campo electromagnético, corresponden a la 
familia de discos de Morgan-Morgan [63, 62]. Con base en estas soluciones, usando métodos estándar 
de generación de soluciones, se pueden obtener modelos de discos con campo electromagnético [36]. 

Una clase más general de soluciones "naturales" que representan discos delgados se obtiene me- 
diante la solución de la ecuación de Laplace (2.3a) en coordenadas esferoidales oblatas con apropia- 
das condiciones de frontera. La solución obtenida para *??(r, z) se usa entonces para obtener la función 
analítica !F(£), necesaria para hacer la transformación (2.6). El problema de valor en la frontera para 
un disco delgado se transforma, por consecuencia, en un problema de valor en la frontera para un 
cascarón en las coordenadas canónicas de Weyl. Así, las soluciones resultantes se pueden interpretar 
como discos delgados finitos o infinitos con presión radial y campo electromagnético. Discos delga- 
dos sin campo electromagnético, obtenidos mediante este mismo procedimiento pueden verse en las 
referencias [62, 13, 28, 22, 69]. 

El presente capítulo expone un formalismo suficientemente general mediante el cual se permite in- 
terpretar ciertas soluciones de las ecuaciones de Einstein-Maxwell en términos de modelos relativistas 
de discos delgados. Con el fin de obtener tal interpretación, en primera instancia, se considerarán solu- 
ciones axialmente simétricas de las ecuaciones de Einstein-Maxwell en las cuales la primera derivada 
del tensor métrico y la primera derivada del potencial electromagnético presenten una discontinui- 
dad finita a través de un disco delgado. Dada entonces una solución con las propiedades anteriores, 
se obtendrán el tensor de energía-momentum superficial y la densidad de corriente superficial de la 
distribución de materia, por medio del formalismo necesario para el tratamiento de las ecuaciones de 
Einstein-Maxwell en términos de distribuciones tensoriales. 

Las propiedades físicas de la distribución de materia se estudiarán entonces mediante un análisis 
detallado del tensor de energía-momentum superficial, considerando soluciones suficientemente ge- 
nerales de las ecuaciones de Einstein-Maxwell que describan una distribución superficial de materia 
no sólo anisótropa, sino que además presente flujo de calor. Posteriormente se presentará una inter- 
pretación alternativa de tales discos con base en el CRM, en el cual el tensor de energía-momentum 
superficial y la densidad superficial de comente se expresan en términos de la superposición de dos 
fluidos perfectos contrarrotantes. 
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2.2. Ecuaciones de Einstein-Maxwell para Cascarones de Mate- 
ria 

Con base en las ideas expuestas se formularán las ecuaciones de campo de Einstein-Maxwell para 
el caso en que existe un cascarón de materia en el espaciotiempo M, de tal forma que las primeras 
derivadas del tensor métrico y del potencial electromagnético presentan discontinuidad finita a través 
de una hipersuperficie S. El tensor métrico g ab y el potencial electromagnético A„ se suponen conti- 
nuos a través de E. Las ecuaciones de campo de Einstein-Maxwell, en unidades geometrizadas tales 
que c = SnG = fio = £o = 1 » están dadas por 

G ab = T llb , (2.8a) 
ud> - i« (2.8b) 



\gabFcdF c \ (2.9a) 
Fot = ¿i», - A,b, (2.9b) 

donde todos los símbolos tienen el significado usual. 

La continuidad en el tensor métrico y en el potencial electromagnético a través de S se puede 
expresar como, 

[fe] = 9 + J* ~ y~a,k = 0, (2.10a) 
[A a ] = AX ~ A;\, = 0, (2.10b) 

mientras que para la discontinuidad en la derivada del tensor métrico se escribe 

y =[&(*)]. (2.H) 

donde ()' denota derivada respecto de (f>. Ahora, usando la aproximación distribucional [68, 58, 87] o 
por las condiciones de junturas de cascarones delgadaos [40, 41, 70] se tiene 

(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 

donde 

' (2.13) 

son los símbolos de Christoffel y se supone que el tensor métrico, g ub , se puede expresar en términos 
de una función de distribución según, 

g al , = (g ab f = cj + ab &(d>) + g-J 1 - &(<p)}, (2. 14) 

siendo &((f>) la función de Heaviside y ó((f>) la distribución de Dirac con soporte en E. 



9ab ~ 


(Qabf, 


9ab,c = 


(g a b,c) D , 


Ifa = 


(r'Lf, 


1 bc,d ~ 




1 be ~ 
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Con el uso de los resultados anteriores en la definición del tensor de curvatura de Riemann, 

«w = r 'L, - Kcm + r M r » - ry*. (2-i5) 

y la suposición que el tensor métrico g a ¡, es, por lo menos, de clase C 3 en las regiones M*, se obtiene 
para el tensor de Riemann la expresión [57] 

«V- = ( R V<,) D + (2-16) 

donde 

«tai = ffU". " FU*. (2-17) 
los (R"^)* son los tensores de Riemann usuales definidos en M* y n a es un vector normal a la 
hipersuperficie S. 

Las discontinuidades en las primeras derivadas del tensor métrico están caracterizadas por el ten- 
sor bab, definido a través de las relaciones 



[gabA = b ab n c , (2.18) 



de tal forma que 



A a 


= (A„) D , 


A aJ , 


= (A aJ ,f, 


pab 




r-ab 


= (F" b b f 



PU = \0>l^ + b%n i -g aí b hc n e \ (2.19a) 
H"tcd = ^(byi b n c -b a c n b n i + b bc n a n i -b ii ri'n c \ (2.19b) 

donde todas las cantidades se evalúan en la hipersuperficie S. 

De igual modo, haciendo uso de la continuidad en el potencial electromagnético A a , se puede 
escribir 

(2.20a) 
(2.20b) 
(2.20c) 

" ,6 H<5(0). (2.20d) 

Suponiendo que el tensor de energía-momentum T ab y la densidad de corriente J" se pueden 
expresar en la forma 

Tab = (T ílb ) D + Qab 5(<Pl (2.21a) 

Ja _ ( ja )D + ¡a (2.21b) 

donde Q a ¡, y I" = ■sj—gl" son el tensor de energía-momentum y la densidad de corriente asociados 
con la hipersuperficie S y los T*. y J" ± , con J" = ■\¡ = g J", son los tensores de energía-momentum y la 
densidad de corriente usuales definidos en M*; se puede probar fácilmente, mediante (2.21b), que las 
ecuaciones de Einstein-Maxwell (2.8) son equivalentes al sistema de ecuaciones 

G% = r¿, (2.22a) 

Hab ~ \gabH = Qab, (2.22b) 

Pfj, = Jl, (2.22c) 
[F ab ]n b = I", (2.22d) 
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donde H ab = H l mb y H = gfH ab . 

Cuando el cascarón de materia es la única fuente de los campos gravitacional y electromagnético, 
de tal forma que el resto del espaciotiempo es vacío (o electrovacío), J" = 0 y las ecuaciones (2.22a) 
y (2.22c) se reducen a las ecuaciones de Einstein-Maxwell en el vacío 

Gf, b = Tf lb , (2.23a) 
Ff. b = 0. (2.23b) 

De esta forma, después de resolver el sistema anterior, las ecuaciones (2.22b) y (2.22d) se pueden 
usar para obtener las expresiones del tensor de energía-momentum y la densidad de corriente para un 
cascarón de materia, 

2£>6 = (K n b ~ bln c )n c + (b c b n c - b c c n b )n" - 6 a b (b e c n e - b e e n c )n c , (2.24a) 
I a = [F" b ]n b , (2.24b) 

donde todas las cantidades son evaluadas en la hipersuperficie X. 



2.3. Discos Delgados Axialmente Simétricos 

En la sección anterior se presentó el formalismo necesario para interpretar soluciones estaciona- 
rias axialmente simétricas de las ecuaciones de Einstein-Maxwell en términos de una hipersuperficie 
X (cascarón) que divide el espaciotiempo M en dos regiones. En la presente sección se restringe 
tal formalismo al caso particular en que la hipersupercie £ a analizar constituye un disco delgado 
axialmente simétrico. Para ello se introduce en el espaciotiempo M el elemento de linea de Weyl- 
Lewis-Papapetrou en coordenadas cuasicilindricas x" = (t, tp, r, z), 

ds 2 = - e m {dt + 'Wdipf + e- m [K 2 d V 2 + e 2A (dr 2 + dz 2 )}, (2.25) 

donde O, A, "W y K son funciones de r y z únicamente y se considera la hipersuperficie £ definida 
por la función (f>(x") = z, con vector normal n„ = (f>,„ = S~. Dada una solución de la forma (2.25) para 
las ecuaciones de Einstein-Maxwell en el vacío el tensor métrico g* b , definido para z > 0, se obtiene 
de la relación 

=g + ab dx a dx", (2.26) 
y el tensor métrico g~ b , definido para z < 0, se obtiene mediante la relación 

íü('-,z) = g*Jr,-z), (2-27) 

de tal forma que [g ab ] = 0. Esto significa construir modelos con simetría de reflexión con respecto al 
plano z = 0. Las relaciones anteriores implican que, cuando z £ 0, 

= -9 + ab ¿r,-z), (2.28) 

así, las discontinuidades en las primeras derivadas del tensor métrico se pueden escribir como 

ba, = [<tej = 2 ./,;,. ! .. (2.29) 
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Para la métrica (2.25), las componentes no nulas del tensor ¿ 0 ¡, son 



b„ = -4e 2 °<D,z, (2.30a) 

b lv = -2e 2 °('W, z + 2'W®-), (2.30b) 

b„ = Ae-™{K(K z -1l<& A )-'W('W- + 'W<í>, z )e 4 ' s '}, (2.30c) 

brr = 4e 2(A - í,) (A, ; - <b.), (2.30d) 

b zz = 4e 2(A - <t) (A, ; - O,-), (2.30e) 



donde todas las cantidades se evalúan en la hipersuperficie z = 0 + . 

Usando los resultados anteriores en la expresión (2.24a), se obtiene para el tensor de energía- 
momentum del disco las componentes no nulas 



Ql ■ 


= ^, {27? 2 (A,--23>,-) + 2KH, Z - e 4 ®1V*Wj , 


(2.31a) 


fií = 


= \2<RrW{K, z - 2ft<b, z ) - {H 2 + WV")OV,.- } , 


(2.31b) 


Ql ■ 




(2.31c) 


Q\ 

Ql ■ 


p2(0-A) 

= ^ 2 (2ft 2 A, z + e 44, 'H"H>«), 

2(«-A) 


(2.3 Id) 
(2.31e) 



aquí, como antes, todas las cantidades se evalúan en la hipersuperficie z = 0 + . 

De otro lado, el "verdadero" tensor de energía-momentum superficial (SEMT) del disco, S„¡,, y la 
"verdadera" densidad de corriente superficial, j,„ se pueden obtener a través de las relaciones 

S ab = f Qab 5(z) ds„ = e"-" Q ah , (2.32a) 

ja = f h 6{z) ds„ = e A -' l I a , (2.32b) 

donde ds„ = sjg^ dz es la medida física de longitud en la dirección normal al disco. 
De modo que, las componentes no nulas de SJ¡ vienen dadas por 



S°o ■ 


O-A 

= -— - (27? 2 (A,- - 2Í>,, ) + 2«,-c 4 'W, ! ), 


(2.33a) 


S° = 


O-A 

= -— - [29rW(%, z - 2ftO, ; ) - (K 2 + r W 2 e 4,t ') r W, z } , 


(2.33b) 


si ■ 


-O-A 

= V^)' 


(2.33c) 


s\ ■ 


= ^-(27? 2 A, ; +e 4<I> nynV, ; ), 


(2.33d) 


si ■ 


O-A 

= ^(2^)- 


(2.33e) 
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y la densidad de corriente en la superficie del disco mediante 



.1, 

h 



(2.34a) 
(2.34b) 



donde se ha hecho uso de la simetría de reflexión los potenciales métricos y electromagnéticos con 
respecto a la superficie del disco, se ha elegido el cuadripotencial electromagnético en la forma A„ = 
(A,, A,/,, 0, 0) y, como antes, todas las cantidades son evaluadas en z = 0 + . 



2.4. Interpretación Física de los Modelos de Discos Delgados 
2.4.1. El Problema de los Valores Propios para el Tensor Energía-Momentum 

Ahora, con el fin de analizar las características del SEMT y las densidades de carga y corriente 
eléctrica, es conveniente expresar éstas en términos de una tétrada ortonormal. Se usará la tétrada del 
"Observador Localmente Estático" (LSO) [43], u observador en reposo con respecto al infinito, el 
cual está dado por 



<D-A iríi 



de modo que 



donde 



„e" e" + S, lim e" e" + 5,„„e° e" - S„ 

» (0) (0) (1KI) (1) (1) (2X2) (2) (2) (1 



e" e" + e" e 



-*} 

(i) (0)' 



2e" 



2<D,-- A,. 



<R 

son las componente no nulas del SEMT con respecto al LSO y, de igual manera, 



/(O, = -2é-- a A,, ; , 



= 2- 



n 



<'W'A- A.\ 



(2.35a) 

(2.35b) 
(2.35c) 
(2.35d) 

(2.36) 

(2.37a) 

(2.37b) 
(2.37c) 
(2.37d) 

(2.38a) 
(2.38b) 



son las densidades de carga eléctrica y corriente eléctrica azimutal del disco medidas por este obser- 
vador. 
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lab 



El caso más simple ocurre si S mm = 0 tal que el SEMT es diagonal. Este caso corresponde 
a una solución estática, "W = 0, en consecuencia se tienen discos delgados estáticos con campo 
electromagnético. También se puede tener S ím) = 0 para soluciones estacionarias con ('1i / , z )| 0+ = 0 y 
entonces esta clase de soluciones se puede considerar como una representación de "discos localmente 
estáticos". En ambos casos se tiene el SEMT de un fluido general, 

eV a V b + Pl X a X b + Pl Y a Y b , (2.39) 

con densidad de energía e = S mm y esfuerzos principales p, = S (¡m y p 2 = 5 (2)(2) , medidos por el 
LSO. La velocidad del fluido viene dada por V" = e" 0) y las direcciones principales son X a = e" y 

ya _ a 

(2) 

De otra parte, cuando S m¡) + 0 el contenido físico de SEMT se puede analizar escribiendo éste en 
la forma canónica. Para tal fin se resuelve el problema de valores propios para el SEMT y se expresan 
las cantidades físicas que le caracterizan en términos de valores y vectores propios. 

A continuación se soluciona el problema de valores propios para el SEMT S ab , para el caso en 
que S m » * 0. Expresando el problema de valores propios en la tétrada ortonormal del LSO, 

S*n?° =A«7 WW Í», (2.40) 

se pueden obtener fácilmente las soluciones Á ± = (B ± VD)/2, á 2 = S 0){2) y A 3 = 0, donde 

+ 5 (0)(0) , (2.41a) 
-S<„xo,. (2.41b) 

(2.41c) 

4C 2 . (2.4 Id) 

Los correspondientes vectores propios están dados por 

(A + y¡D\ ) 

(2.42b) 
(2.42c) 

donde los N ± son factores de normalización apropiados. 

Nótese que existen tres diferentes casos dependiendo del valor de D. En el primer caso, cuando 
D > 0, los dos valores propios A ± son reales e iguales y los dos vectores propios son reales 
y ortogonales, así uno es espacial y el otro temporal. Escogiendo apropiadamente los factores de 
normalización N ± , es fácil ver que 

f ± -& = ±v,, (2.43) 

donde v A = A/\A\, de modo que v A determina el carácter causal de los vectores propios f ± y se obtiene 
una tétrada ortonormal con ellos y los otros dos vectores. En el segundo caso, cuando D < 0, los dos 
vectores propios f ± son complejos conjugados y ortogonales. Se pueden escoger, apropiadamente, los 
factores de normalización N ± de tal manera que 

í + -£ + =f_-í_<0, (2.44) 



A 


_ c 

u mm 


B 


- 

u mm 


C 


- c 

u w>m 


D 


= A 2 - 



ñ 


= N ± 


e° - 

I W) 




= Y" 


= e" 

(2) 


ñ 


= Z" 


= e" 

(3) 
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y entonces los vectores propios se pueden escribir como 

=£ + ¿ff, (2-45) 

con 

&•& < 0, (2.46a) 

> 0, (2.46b) 

= 0. (2.46c) 

Normalizando las partes real e imaginaria de los vectores propios, se tiene nuevamente una tétrada 
ortonormal formada con t¡ R , £;/ y los otros dos vectores propios. Finalmente, D = 0, A + = A_ y los 
dos vectores propios £ ± son colineales y nulos, ¿¡. ± ■ f ± = 0. Estos vectores nulos se pueden expresar 
apropiadamente en términos de un vector temporal y uno espacial y se tiene nuevamente una tétrada 
ortonormal. Además, se puede probar fácilmente que todas las cantidades que caracterizan el tensor 
energía-momentum son funciones continuas de D. De acuerdo con ello, se puede considerar el caso 
D = 0 como el límite de los casos anteriores cuando D — > 0. 

En resumidas cuentas, en los tres casos se obtiene una tétrada ortonormal asociada con los vectores 
propios. Esta tétrada se puede escribir apropiadamente como 



donde 



V a 








Vi-ü 2 


X" 


= -v 








i Vi - V 2 


ya 


= e" , 

(2)' 




z a 


= e O)' 





( (v„ VI>-A)/2C; L>>0, 
\ (v c V 3 » - 2C)/A; D < 0, 



(2.47a) 

(2.47b) 

(2.47c) 
(2.47d) 

(2.48) 



con v c = C/\C\, v = 1 cuando D > 0 y v = v c cuando D < 0. 

La tétrada ortonormal [V , X" , Y" , Z"} es comóvil con el disco, tal que el vector velocidad del 
disco está dado por el vector propio temporal V". De acuerdo con ello, la velocidad "tangencial" o 
"circular" del disco v, medida por el LSO, está dado por (2.48). Ahora, es fácil ver que en todos los 
casos se tiene que 

\v\ < 1, (2.49) 

y que 

lím|«|=l. (2.50) 

D->0 

Esto es, la velocidad tangencial de los discos es siempre menor que la de la luz y solo alcanza este 
valor cuando D = 0, el caso cuando se tiene un vector propio nulo. 
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2.4.2. Contenido Físico del Tensor Energía-Momentum 

Expresado en términos de la tétrada comóvil {V a , X a , Y", Z"), el SEMT se puede escribir en la 
forma canónica 

Sal, = eVaVi+p.XaXt+pJaYt + qfVaXt+XaV,,). (2.51) 
Entonces el SEMT representa un fluido general con flujo de calor diferente de cero. La densidad de 
energía del disco está dada por 

(v A Vo-fi)/2, D>0, 

-B/2, D < 0, K ' 
mientras que 

r (n VD + B)A D>0, 

' 1 \ B/2, D < 0, y ' 
es el esfuerzo principal en la dirección de X a , y 

Pi = S m) (2.54) 

es el esfuerzo principal en la dirección de Y". El vector flujo de energía está dado por 5" = —S a y h . 
Entonces, usando la expresión (4.4) para el SEMT y la ortonormalidad de la tétrada comóvil, se puede 
escribir 

S" = eV" + qX", (2.55) 
de modo que el flujo de energía tiene dos contribuciones: la del "flujo de energía cinética" 

e" = eV, (2.56) 

y la del vector flujo de calor, 

q" = qX", (2.57) 

donde 

0, D > 0, 

V^D/2, D < 0, (2 ' 58) 
es la función de calor. 

La condición débil de energía requiere que la densidad de energía sea no negativa, 

€ > 0, (2.59) 

mientras que la condición de energía fuerte es equivalente a la condición 

p > 0, (2.60) 

donde 

p = € + P¡ +p 2 (2.61) 
es la "densidad newtoniana efectiva". Finalmente, la condición dominante de energía requiere que el 
vector flujo de energía sea temporal o nulo, 

S"S„ < 0. (2.62) 
El SEMT (4.4) satisfará las condiciones anteriores cuando 

\q\ < leí, (2.63) 
garantizándose así la propagación causal de la energía. 
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2.5. El Modelo de Contrarrotación 

Ahora se considerará, con base en las referencias [20], [54] y [16], la posibilidad de que el SEMT 
S ab y la densidad de corriente j" se puedan expresar como la superposición de dos fluidos cargados 
contrarrotantes que circulan en direcciones opuestas. Así, se supone que S y j" se pueden escribir 
como 

S ab = Sf + S'L\ (2.64a) 
/ = /I + £, (2.64b) 

con el SEMT y la densidad superficial de corriente de los fluidos contrarrotantes prógrados y retró- 
grados dados, respectivamente, por 

Sf = (e + + p + ) U a + U\ + p + h" b , (2.65a) 

Sf = (c + p_) U"_U b _ + p. h" b , (2.65b) 

j+ y j-, donde 

ft°* = g a *-e;e* )) (2.66) 

es la métrica inducida sobre la hipersuperficie z = 0. 

Con el fin de descomponer (2.64a), inicialmente se escribe el SEMT como 

S" b = {S mm + S_J e" e b + (5 ...... - S_J e" e* - 5, 0H „(e" e 6 + e" e* } + S m2 ,h ab . (2.67) 

> (OK0) (2)(2)J (0) (0) 1 (D(l) (2)(2)J (! ) (i) (0)(1) l (0 ) (l) (l) (0) ' (2)(2) V / 

Ahora, si = (í/°, U±, 0, 0) son respectivamente los vectores velocidad de los dos fluidos contrarro- 
tantes, su proyección sobre la tétrada del LSO da 

e" + u+s" 

m = ' , (2.68a) 

Vi - «+ 

e" + U-t" 

U" = '"' " ' . (2.68b) 
Vi - u 2 _ 

donde u ± = U+' /1¡T son las velocidades tangenciales de los fluidos contrarrotantes con respecto a la 
tétrada LSO. 

Ahora, expresando los vectores de la tétrada e° y e° en términos de los vectores velocidad U" y 
[/f mediante el uso de (2.68a) and (2.68b), se obtiene 

S<A = f(u.,u.)(l-u 2 + ) uaub + f(u + ,u + )(l-u 2 _) u<¡ub 
(u + - u_) 2 + + (w + - í¿_) 2 

f(u + , »_)[(! - ul)(l - u 2 _)]' 



(u + - íí_) 2 



U a + Ut + U a _ U° + + S ,,„,,/!"", (2.69) 



donde 



f(ui, u-i) — [S (0)(o) + S mf2) }uiii2 + S m¡t [u¡ + 112} + {S ím ¡ ~ S & d 

}. (2.70) 

entonces, con el fin de escribir el SEMT en la forma (2.64a), los términos mixtos tienen que anularse 
y las velocidades tangenciales de contrarrotación estarán relacionadas por 

/(«+,«_) = 0, (2.71) 
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donde se asume que u± ¿ 1. Así, el SEMT se puede escribir como la superposición de dos fluidos 
cargados si, y solo si, la ecuación anterior admite una solución tal que u + + «_ . Este último resultado 
generaliza la obtenida en [22] y es completamente equivalente a la condición necesaria y suficiente 
obtenida en [16]. Sin embargo, la relación anterior no determina u ± únicamente y por consecuencia 
existe cierta libertad en la escogencia de los vectores U%. 

De otra parte, asumiendo un escogencia dada para las velocidades de contrarrotación en acuerdo 
con la ligadura (2.71), se puede escribir el SEMT y la densidad de carga eléctrica de los fluidos como 
(2.64a) con 



e+ + P+ 

e_ + p_ 
P+ + P- 
o-± 



1-KÍ 



U- — U+ 
l-u 2 



U+ — U- 

W (2){2)> 



[(■^(0X0) + ^(2)(2)í"+ + S<! 



"±/-,<lj _ (<»„ 



u± — 



y, con las expresiones anteriores, adicionalmente se obtiene la relación 



+ S m 



(2.72a) 

(2.72b) 
(2.72c) 

(2.72d) 
(2.73) 



Estas expresiones son la generalización estacionaria de la obtenida en [22] para el caso estático y 
también son equivalentes a las expresiones en [16]. Nótese que la densidad de energía, e ± , y la presión, 
p ± , no están únicamente definidas por las relaciones anteriores, tampoco para valores definidos de u ± . 
Acordemente, se necesita imponer alguna condición adicional sobre la densidad de energía y presión 
con el fin de determinar completamente el CRM. 



2.5.1. Determinación de las Velocidades de Contrarrotación 

Tal como se anotó en la sección anterior, es necesario suponer alguna relación adicional entre las 
velocidades de contrarrotación con el fin de determinar completamente el CRM. La más simple, y 
comúnmente considerada, posibilidad es tomar las dos velocidades tangenciales de contrarrotación 
como iguales y opuestas; esto es, 

u+ = — u- = u, (2.74) 
donde u es obtenida a partir de (2.71) y viene dada por 

2 _ jVjOj ~ 3(2X2) _ KA.; ~ ft¿ 

Ahora, con el fin de que las expresiones anteriores tengan sentido, se necesita imponer la condición 
adicional 

0<m 2 <1, (2.76) 

ya que, si esta última condición no es satisfecha, no se puede obtener un CRM bien definido. Sin 
embargo, puede pasar que la solución considerada satisfaga las condiciones anteriores sólo en una 
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región parcial del disco y, si este es el caso, la correspondiente interpretación del CRM es válida sólo 
parcialmente. 

Otra posibilidad comúnmente aceptada es suponer que los dos fluidos se mueven a lo largo de 
electrogeodésicas. Con el fin de analizar esta posibilidad, se escriben los vectores velocidad como 
U% = t/!J(l, Í2 ± , 0, 0), de tal modo que la ecuación de la electrogeodésica es (para ver detalles de la 
obtención de esta ecuación ver Apéndice A) 

\(U*) 2 (g m . r + 2g 0l ,Q. ± + g n M)(e± + P±) ~ P±,r = -°-±U° ± + <V n *) > ( 2 - 77 ) 

donde, como se sabe, Í2+ = U+/U' ± son las velocidades angulares de los fluidos contrarrotantes. 

En conclusión, las velocidades tangenciales de contrarrotación pueden ser determinadas explíci- 
tamente si y sólo si se asume una condición adicional entre ellas, como por ejemplo que estas sean 
iguales y opuestas o la condición de movimiento de los fluidos a lo largo de electrogeodésicas. Aho- 
ra, puede pasar que las soluciones obtenidas no satisfagan ninguna de estas condiciones. Esto es, las 
velocidades de contrarrotación, en general, no son determinadas completamente mediante las ligadu- 
ras (2.71). Entonces, el CRM es en general indeterminado dado que las densidades de energía y las 
presiones no pueden ser escritas explícitamente sin un conocimiento previo de las velocidades tan- 
genciales de contrrrotación. En efecto, como se muestra en [16] y [37], se puede interpretar el SEMT 
en términos de una familia 2-paramétrica de fluidos perfectos. 



2.5.2. Determinación de las Presiones y Densidades de Energía de Contrarro- 
tación 

Supóngase ahora que las velocidades tangenciales de contrarrotación se conocen a través de una 
de las dos vías (o mediante otro camino) descritas en la sección anterior, se obtendrán entonces las 
expresiones para las densidades de energía y las presiones de los dos fluidos contrarrotantes. Como 
ya se dijo, el sistema (2.72a) - (2.72c) es, en general, indeterminado. Así que no es posible obtener 
las densidades y presiones sin una condición adicional. No obstante, si S (2)(2) = 0, se puede obtener 
una solución del sistema (2.72a) - (2.72c) bajo la condición que los "esfuerzos" de contrarrotación 
sean positivas (presiones). Así, se tiene que las presiones son 

P + = P- = 0, (2.78) 

de tal modo que el SEMT (2.36) puede interpretarse como el correspondiente a dos fluidos de polvo 
contrarrotantes. 

Con los resultados anteriores, las expresiones (2.72a) y (2.72b) se reducen a 

e + = S m« u - + S ™ (2 . 79a) 
(m_ - u + )(l — u%)~ 



S (lm U + + 5(0X1) 

(u + - »_)(! - u 2 _)~ l ' 



(2.79b) 



y entonces todas las cantidades que caracterizan el CRM están completamente determinadas. Como 
se puede ver, las densidades de energía de los dos fluidos de polvo contrarrotantes son diferentes, 
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también cuando se tiene iguales y opuestas velocidades tangenciales. Sin embargo, si S mn = 0 y 
las velocidades tangenciales son iguales y opuestas, las densidades de los dos fluidos de polvo son 
iguales y opuestas. 

De otro lado, cuando 5 121(2I £ 0, se tiene que imponer una condición adicional con el fin de obtener 
explícitamente las expresiones para las densidades y presiones. Una posible, y común condición es 
suponer que las densidades de energía de los dos fluidos contrarrotantes son iguales. Así, la ecuación 
(2.73) conduce a 

^ = ^ = ^S mm -S mm +S aia }, (2.80) 
mientras que las ecuaciones (2.72a) - (2.72b) conducen a 

P+ = { s ,mi> ~ [S<oxo) + 5 o»)] K«+ + «-)/("+ - «-)]} /2 - S mm [(1 + ¡í+M_)/(tí + - «_)] , (2.81a) 

p_ = {S im + [S m) + S m) ] [(u + + «_)/(»+ - «_)]} 12 + S m) [(1 + « + «_)/(« + - «_)] . (2.81b) 

Así, las presiones de los dos fluidos contrarrotantes son diferentes, también cuando las velocidades 
tangenciales son iguales y opuestas. Las presiones son iguales sólo si S, m , = 0 y las velocidades 
tangenciales son iguales y opuestas. 

Ahora, en lugar de suponer que las dos densidades son iguales, se puede suponer que las presiones 
son iguales, así que la ecuación (2.72c) conduce a 

p + = P- = S atm /2, (2.82) 

y las ecuaciones (2.72a) - (2.72b) conducen a 

e + = [S w / 2] -[($ im u^ + S m) u + )/(u + -U-)]-S mi [(l + u + uJ/(u + -»_)], (2.83a) 

e_ = [5 (2>21 /2] + [(5 (()K0) ¡ í _ + 5 ( , K1 ,» + )/( M+ -í<_)]+5 mil ,[(l + í< + i«_)/(« + -t(-)]. (2.83b) 

Así, se tiene entonces dos fluidos contrarrotantes con diferentes densidades de energía, también cuan- 
do las velocidades son iguales y opuestas y, como en los previos casos, las densidades de energía son 
iguales sólo si 5 (0KI1 = 0 y las velocidades tangenciales son iguales y opuestas. 
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3.1. Introducción 

continuación se presenta una familia infinita de discos de polvo cargados eléctricamente. Estos 
discos, axialmente simétricos de extensión finita, se obtienen solucionando las ecuaciones de 
Einstein-Maxwell para espaciotiempos conformestáticos, los cuales se caracterizan por sólo una fun- 
ción métrica. Con el fin de obtener las soluciones, se establece una relación funcional entre la función 
métrica y el potencial eléctrico. También se supone que la función métrica depende de otra función 
auxilar, la cual se toma como una solución de la ecuación de Laplace [33]. Las soluciones para la 
función auxiliar elegidas corresponden a las dadas por la familia de discos generalizados de Kalnajs 
obtenidas por González y Reina [30], expresadas en términos de las coordenadas esferoidales oblatas, 
las cuales corresponden a una familia de discos newtonianos delgados axialmente simétricos de radio 
finito. Los discos relativistas delgados obtenidos tienen densidad de carga que es igual, excepto por un 
signo, a su densidad de masa. De modo que la densidad de energía de los discos es siempre positiva y 
bien comportadas, anulándose en el borde del disco. Acordemente, como los discos obtenidos están 
conformados por polvo, su tensor energía-momentum está de acuerdo con todas las condiciones de 
energía. El presente capítulo corresponde a la versión en lengua española del trabajo [27]. 

3.2. Ecuaciones de Einstein-Maxwell y Soluciones tipo Discos Del- 
gados Finitos Conformestáticos 

Como es sabido, para la descripción de un espaciotiempo estacionario axialmente simétrico puede 
tomarse el elemento de linea de Weyl-Lewis-Papapetrou en coordenadas canónicas de Weyl dado por 



Capítulo 3. Discos Finitos de Polvo, Eléctricamente Cargados, en Espaciotiempos 

Conformestáticos Axialmente Simétricos 

(2.5). Ahora bien, si se hace la transformación 

ft^r, (3.1a) 
Z^z, (3.1b) 

y en (2.5) se escoge IT = 0, para la función métrica ü, mientras que para la otra función métrica 
¥ se toma ¥ = A, donde A es función sólo de r y z, se tiene el tensor métrico correspondiente a 
un espaciotiempo comformestacionario axialmente simétrico. En el caso particular en que JK = 0 se 
obtiene entonces el elemento de linea 

ds 2 = - e 2Á dt 2 + e- 2 \r 2 díp 2 + dr 2 + dz 2 ), (3.2) 

el cual corresponde a un espaciotiempo conformestático axialmente simétrico. Este tipo de métricas 
pueden interpretarse geométricamente como métricas en las cuales la parte espacial es conformalmen- 
te plana o, dicho en otros términos, que su parte espacial es obtenida a partir de una transformación 
conforme de la parte espacial de la métrica minkowskiana. Como se verá más adelante, este tipo de 
métricas corresponden al exterior de configuraciones constituidas por polvos estáticos eléctricamente 
cargados en los cuales la fuerza eléctrica y la gravitacional se encuentran en perfecto balance [86]. 

Ahora bien, para la métrica 3,2 el sistema de ecuaciones de Einstein-Maxwell en el vacío en 
unidades geometrizadas tal que c = G = /¿o = £o = 1 es equivalente a 

W 2 A = e- 2Á V4>-V(f>, (3.3) 
V 2 0 = 2VA-V(¡>, (3.4) 

A.,A.j = e^fo </>j, (3.5) 

donde los índices i, j = 1 , 2, 3 y V es el operador diferencial usual en coordenadas cilindricas y se 
ha escogido el potencial electromagnético en la forma 

A, = (A, = -<f>, 0, 0, 0), (3.6) 

y se ha supuesto que tanto la función métrica como el potencial eléctrico son independientes de í. 

Con el propósito de solucionar el sistema de ecuaciones anterior, primero se esupone que el poten- 
cial eléctrico ¡p es funcionalmente dependiente de la función métrica A, (f> = (f>(A), tal que la ecuación 
(3.5) implica que 

[<P'(A)] 2 = e 2 \ (3.7) 

cuya solución está dada por 

(f> = ±e Á + k,, (3.8) 

con k¡ una constante de integración arbitraria. Bajo esta solución, el sistema de ecuaciones (3.3)-(3.4) 
se reduce únicamente a la ecuación diferencial parcial no lineal, 

W 2 A = VA ■ VA, (3.9) 
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para la función métrica A. Si ahora se asume una dependencia funcional para la función métrica A{U), 
donde U es una función auxiliar elegida como una solución de la ecuación de Laplace, la ecuación 
(3.9) se reduce a 

A"(U) = [A'(U)]\ (3.10) 

cuya solución viene dada por 

^ = 77^. (3.1D 
U + k 2 

donde k^ y £3 son constantes de integración. 

Como es natural, si se espera un comportamiento físicamente apropiado en el infinito, se imponen 
algunas condiciones en las soluciones. Así, con el fin de obtener un espaciotiempo asintóticamente 
plano, se exige que e Á = 1 en el infinito. También es necesario que 0 = 0 en el infinito, dado que 
se pretende obtener fuentes de extensión finita. Consecuentemente, se considerarán para la función 
U sólo aquellas soluciones de la ecuación de Laplace que exhiban el comportamiento U = 0 en el 
infinito. De esta manera, tomando k^ = £3 se obtiene el comprtamiento deseado para A. Finalmente, 
para completar las exigencias impuestas sobre el comportamiento de 0, se toma k¡ = +1, así que la 
ecuación (3.8) puede ser escrita como 



0: 



U + k 



(3.12) 



donde las constantes han sido renombradas como k^ = k¡ = k. 

Como se puede ver de (3.11) y (3.12), las soluciones del sistema de ecuaciones de Einstein- 
Maxwell (3.3)-(3.4) son expresadas en términos U, una solución de la ecuación de la ecuación de 
Laplace. De manera que, con el objetivo de tener soluciones correspondientes a discos delgados finitos 
a se considerarán solamente soluciones de la ecuación de Laplace que describan propiamente fuentes 
discoidales Newtonianas de radio finito. Esta clase de soluciones se pueden obtener introduciendo 
las coordenadas esferoidales oblatas, las cuales se adaptan en forma natural a la fuente debido a la 
geometría de estas fuentes. Estas coordenadas están relacionadas con las coordenadas cilindricas a 
través de las expresiones [64] 

r 2 = a 2 (l+f 2 )(l-/, 2 ), (3.13a) 
Z = a%r], (3.13b) 

donde 0<£<ooy-l<;/<l.EI disco tiene coordenadas f = 0, 0 < rf < 1 y, a través del disco, r¡ 
cambia de signo pero no cambia su valor absoluto. 

El comportamiento singular de la coordenada 77 implica que un polinomio de r¡ con potencia par 
siempre es continuo en cualquier parte pero presenta una discontinuidad en las derivadas con respecto 
de r¡ en el disco. Acordemente, la solución general de la ecuación de Laplace correspondiente a una 
fuente tipo disco delgada axialmente simétrica de radio finito a puede expresarse como [2] 

V(g, rj) = - ¿ C lnqi M)PiÁn\ (3.14) 

donde los P2„(v¡) son los polinomios de Legendre de orden 2n y ¡^«(f) = i 2 "* QiJM), con Q2A'?) 
las funciones de Legendre de segunda clase con argumento imaginario. Los C2,, son constantes que 
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habrán de ajustarse apropiadamente de acuerdo a la solución en particular y están directamente re- 
lacionadas con la densidad de masa Newtoniana de la fuente considerada. En la próxima sección 
se presenta una escogencia particular de estas constantes correspondientes a una familia infinita de 
discos con una densidad superficial de masa bien comportada, la familia de discos generalizados de 
Kalnaj recientemente presentados por González y Reina [30]. 

Ahora bien, para la métrica (3.2), la única componente diferente de cero Q' b es 

(3.15) 



S¡ = -^, (3.16) 



2n 

por lo tanto, la única componente diferente de cero del tensor energía-momentum superficial S" b está 
dado por 

e% 
2n 

donde todas las cantidades son evaluadas en z = 0 + . 

El tensor energía-momentum superficial del disco y la densidad superficial de corriente pueden 
escribirse como 

S ab = eV"V h , (3.17a) 



/ = a-V (3.17b) 

donde 

V" = ¿"''(1, 0,0,0) (3.18) 

es el vector velocidad de la distribución de materia. Como consecuencia de lo anterior, la densidad de 
energía y la densidad de carga eléctrica de la distribución de materia están dadas por 

e x X_, 



2n 



2n 

respectivamente. Ahora, mediante el uso de (3.8), la expresión (3.20) puede ser escrita como 



(3.19) 
(3.20) 

(3.21) 



tal que la densidad de carga eléctrica de los discos es igual, excepto por un signo, a su densidad de 
masa. Por consecuencia, las fuerzas eléctrica y gravitacional están en perfecto balance, como en las 
configuraciones de ECD; tal como se mencionó en la introducción al capítulo. 

Ahora, en el contexto de la relatividad general es necesario que el tensor energía-momentum 
satisfaga ciertos condiciones enmarcadas en las condiciones de energía débil, fuerte y dominante 
[39]. En efecto, para el caso de discos de polvo, todas estas condiciones se reducen a la condición que 
la densidad de energía es mayor o igual que cero, e > 0. De otra parte, de la ecuación (3.11), se tiene 
que la densidad de energía puede ser escrita como 
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donde 

U. 

£ = -í (3.23) 

es la densidad de masa Newtoniana de una fuente discoidal con potencial gravitacional dada por U. 
Acordemente, si la densidad de masa Newtoniana X es no negativa en todas partes, la correspondiente 
densidad de energía relativista € será también no negativa en cualquier parte sólo si k < 0. Así enton- 
ces, como se quiere que el tensor enrgía-momentum esté de acuerdo con as condiciones de energía, las 
constantes Ci„ en C2,, deben de escogerse apropiadamente de tal modo que X > 0. Además, si se tiene 
un potencial newtoniano U negativo en cualquier parte, tal como se espera para fuentes Newtonianas 
compactas, entonces la densidad de enegía e será no singular en el disco. 



3.3. Una Familia Particular de Discos 



Ahora se restringirá el modelo general estudiado en la sección anterior considerando una familia 
particular de discos con densidad superficial de energía bien comportada. Los miembros de la familia 
de discos se expresan en términos de soluciones particulares de la ecuación de Laplace U m obteni- 
das escogiendo apropiadamente las constantes de la solución general (3.14). La solución particular 
representa el potencial gravitacional de un disco delgado finito con densidad de masa dada por 



(2m+ 1)M 



2na 2 



r 2 



(3.24) 



donde M y a son, respectivamente, la masa total y el radio del disco y se tiene que m > 1. Para cada 
valor de m, las constantes Cj n están definidas a través de la relación [30] 



donde 



Ki„ 



M 
2a 



(2n+i) ?2n+1 (0)' 
t 1/2 (4n + 1) (2m + 1)! 



(3.25) 



(3.26) 



2 2m (m - n)\T(m + n + |) 

para n < m, y C2,, = 0 para n > m. En la tabla 3.1 se listan los valores diferentes de cero para 
los primeros seis miembros de la familia. Así, usando estas constantes en la solución general (3.14), 
se puede ver que el potencial gravitacional U m será no negativo siempre, tal como se impuso en las 
secciones previas. 

Con los valores anteriores para las contantes C2,,, se puede calcular fácilmente la densidad de 
energía de los discos usando la ecuación (3.22). Ahora, con el fin de ilustrar gráficamente el compor- 
tamiento de los diferentes modelos particulares, primero se introducen cantidades adimensionales a 
través de las relaciones 

aU„,(f) 



U m (f) = 



M 



na 2 I.„,(f) 
M ' 



(3.27) 



(3.28) 



na€ m (r), 



(3.29) 
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Cuadro 3.1: Las constantes C2,, para m = 0, 6. 
m Co C2 C\ Cb C$ Cío Cn 



1 


M 


K 










a 












2 


M 


10AÍ 


3 Ai 








a 


7o 


7a 








3 


jU 


5A4 


9,v; 


5 Ai 






a 


3o 


1 la 


33a 






4 


M 


20M 


162AÍ 


4M 


7M 




a 


Ha 


143a 


1 la 


143a 




5 


a 


25 Ai 
13a 


18AÍ 
13a 


10AÍ 
17o 


35AÍ 
247a 


63AÍ 
4199a 


6 


M 


2M 


27M 


260M 


5M 


378AÍ 


a 




17a 


323o 


19a 


7429a 



33AÍ 
7429o 



donde f = r/a, O < f < 1, y £/„,(?) es evaluada en z = 0 + . 

Consecuentemente, las densidades adimensionales de energía e(P) pueden escribirse como 

e m (r) = - - , (3.30) 

[[/,„(r) + fe] 2 

con A: = (ka)/M. 

Así, usando las anteriores expresiones y los valores de las constantes C2,, dados en la tabla 3.1, se 
obtienen las siguientes expresiones 



3feVT 



(3.31) 
(3.32) 
(3.33) 



£4 = = 5= , (3.34) 



2[fe+ f(? 2 -2)] 2 ' 

5fe(l -r 2 ) 3 ' 2 
2[fe- HO^-S^ + g)] 2 ' 

7it(l - r 2 ) 5 ' 2 
2[k + §§ (5r 6 - 18r> + 24f 2 - 16)] 2 ' 

9fe(l - r 2 ) 7 ' 2 

2[fe - ^|(35r 8 - lóOr 6 + 288F 4 - 256r 2 + 128)] 2 ' 

llfe(l -r 2 ) 9/2 

2 [fe + y|fif2(63í" 0 - 350r 8 + 800r 6 - 960P 1 + 640r 2 - 256)] 2 ' 



(3.35) 



13fe(l-r 2 )"/ 2 

e 6 = = t^tí (3.36) 

2[fe - ¿f^f^lr 12 - 1512r 10 + 4200r 8 - 6400F 6 + 5760? 14 - 3072r 2 + 1024)] 2 

para las densidades de energía de los seis primeros modelos. 
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En la figura 3.1 se gráfica, la densidad de energía superficial adimensional e m como una función 
de f para los primeros seis modelos con m = 1, 2, 3, 4, 5 y 6. En cada caso, se gráfica e m (f) para 
0 < f < 1 con diferentes valores de los parámetros k. Primero se toma k = -0,25, la curva en la parte 
inferior en cada cuadro, y luego k = -0,5, -1, -2, -4 y -8. Como se puede ver, en todos los casos 
la densidad de energía es positiva en todas partes y se anula en el borde del disco. Nótese que exhibe 
diferentes comportamientos dependiendo de los valores de m y k. Así, para los primeros dos modelos 
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Figura 3.2: Densidad de energía superficial e m en función de f con k = -1 para los seis modelos de 
discos con m = 1, ... , 6. 

con m = 1 y m = 2, se encuentra que para pequeños valores de \k\ la densidad de energía presenta 
un máximo cerca el borde del disco, mientras que para valores grandes de \k\ el máximo ocurre en 
el centro del disco. De otra parte, para m > 3 se encuentra que los valores de \k\ el máximo de la 
densidad de energía ocurre en el centro del disco. 

También se puede ver que como los valores m se incrementan, la densidad de energía se concentra 
más en el centro del disco. También, conforme m aumenta el valor de la densidad de energía en la 
parte central de los discos se incrementa de modo que, para un valor dado de k, la densidad de energía 
central es mayor para valores mayores de m. De otra parte, en el borde de los discos el comportamiento 
es opuesto, los valores de la densidad de energía decrecen conforme m se incrementa. Ahora, con la 
intención de ver esto, en la figura 3.2 se gráfica la densidad de energía adimensional ?,„ como una 
función de f con k = - 1 para los modelos de discos con m = 1, ... , 6. Finalmente, para un valor dado 
de m, el valor de la densidad de energía crece con el valor de pero luego decrece conforme|£| se 
incrementa. Esto se describe en la figura 3.3 donde se gráfica la densidad de energía adimensional en el 
centro de los discos, e(0), como función de k para los mismos seis modelos previamente considerados. 
Ahora bien, como la densidad de carga de los discos es igual, excepto por un signo, a su densidad de 
energía, todos los análisis previos son también aplicables al comportamiento de la densidad superficial 
de carga de los modelos. 

A manera de conclusión del tema estudiado en el presente capítulo se puede decir que, se ha 
presentado una familia infinita de discos axialmente simétricos de extensión finita con densidad su- 
perficial de energía y carga eléctrica bien comportados. Los discos fueron obtenidos solucionando el 
sistema de ecuaciones de Einstein-Maxwell en el vacío para un espaciotiempo conformestático. Con 
el fin de obtener las soluciones, se asumió una dependencia funcional entre el potencial eléctrico y la 
función métrica y entre esta y una función auxiliar, las soluciones fueron entonces expresadas en tér- 
minos de una solución de la ecuación de Laplace correspondiente a una familia de discos newtonianos 
delgados de radios finitos, los discos de Kalnajs generalizados [30]. 

Los discos delgados relativistas aquí presentados poseen una densidad de carga eléctrica que es 
igual, excepto por un signo, a su densidad de energía, constituyendo así un ejemplo de los comúnmen- 
te llamadas "discos contrapuestos eléctricamente". La densidad de energía de los discos es siempre 
positiva y bien comportada, cayendo a cero en el borde. También, conforme los valores de m se incre- 
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Figura 3.3: Densidad de energía superficial en el centro de los discos, e(0), en función de k para los 
mismos seis modelos de discos previamente considerados. 

mentan, la densidad de energía es más concentrada en el centro de los discos, con un máximo en r = 0 
para todos los valores de \k\. No obstante, para los dos primeros modelos con m = 1 y m = 2, para 
valores pequeños de \k\ la densidad energía presenta un máximo cerca al borde del disco, mientras 
que para valores grandes de \k\ el máximo ocurre en el centro del disco. 

Además, como la densidad de energía de los discos es siempre positiva y los discos están consti- 
tuidos por polvo, todos los modelos están en completo acuerdo con todas las condiciones de energía, 
un hecho de particular relevancia en el estudio de modelos relativistas de discos delgados. En efecto, 
tal como se mencionó en la introducción al capítulo, muchos de los modelos relativistas de discos 
delgados estudiados en la literatura no está en completo acuerdo con estas condiciones. 

Como se puede ver de las ecuaciones en la sección 3.2, el procedimiento aquí expuesto puede ser 
aplicado no sólo para obtener espciotiempos conformestáticos axisimétricos sino también puede usar- 
se para obtener soluciones no-axilmente simétricas de las ecuaciones de Einstein-Maxwell en el vacío. 
Por consecuencia, los modelos aquí expuestos pueden ser generalizados tomando para la función au- 
xiliar U soluciones de la ecuación de Laplace sin la simetría axial impuesta. Esta idea sigue el curso 
de una trabajo pendiente. Del mismo modo, la generalización a espaciotiempos conformestacionarios 
con campo magnético. 
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n el capítulo a continuación, una familia de nuevas soluciones exactas de las ecuaciones de 



K^J Einstein en el vacío para espaciotiempos estáticos axialmente simétricos es presentada. Toda 
las funciones métricas de las soluciones son calculadas explícitamente y las expresiones obtenidas 
son simples escritas en términos de las coornadas esferoidales oblatas. Además, las soluciones son 
asintóticamente planas y regulares en todas partes, tal como se muestra calculando los escalares de 
curvatura. Estas soluciones describen una familia infinita de discos de polvo delgados con un borde 
interno central, cuyas densidades de energía son siempre positiva y bien comportada, de tal modo que 
su tensor energía-momentum satisface todas las condiciones de energía. Ahora, aunque los discos son 
de extensión infinita, todos tienen masa finita. La superposición del primer miembro de esta familia 
con un agujero negro de Schwararzschild fue presentada en [31]. 

También se presenta una familia infinita de soluciones exactas, completamente integrada explí- 
citamente, de las ecuaciones de Einstein correspondiente a la superposición de un disco de polvo 
contrarrotante con un agujero negro. La solución obtenida representa un disco anular delgado infinito 
(un disco con un borde interno) alrededor de un agujero negro de Scchwarzschild. La masa del disco 
es finita y el tensor de energía-momentum cumple con todas las condiciones de energía. Además, la 
masa total del disco en presencia del agujero negro es menor que cuando el disco está solo. La solu- 
ción puede interpretarse también como si describiera un disco delgado constituido por dos fluidos de 
polvo contrarrotantes que satisfacen también las condiciones de energía. 

Con el fin de describir un espaciotiempo estático axialmente simétrico con un disco infinitesimal- 
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mente delgado como fuente, se hace la transformación 







r, 


Z(r,z) 




z, 






0, 


¥(.r,z) 




$(/, z), 


n(r,z) 




A(r,z), 



que deja el elemento de Weyl-Lewys-Papapetrou (2.5) en la forma [84] 

ds 1 = - e m dt 1 + íT 2< VíV 2 + e 2A (dr + dz 2 )], (4.1) 

en donde las coordenadas x" = (í, <p, r, z) son las usuales coordenadas cilindricas, y se hace uso del 
formalismo general descrito en el Capítulo 2 según el cual existe en el espaciotiempo un disco in- 
finitesimalmente delgado, localizado en la hipersuperficie z = 0, de modo que las componentes del 
tensor métrico son funciones simétricas de z y que sus primeras derivadas con respecto a z poseen una 
discontinuidad finita en z = 0. 

Para esta métrica, las únicas componentes no nulas del tensor energía-momentum Q„t, son 

2e 2 («'- A '[A^-20 )2 ] = Q' n (4.2a) 

donde todas las cantidades son evaluadas en z = 0 + . Ahora bien, usando la tétrada ortonormal 

V" = e-°Sf¡, (4.3a) 

W" = e®<5J/r, (4.3b) 

X" = e°" A ^, (4.3c) 

Y" = e^ó'^ (4.3d) 

se puede escribir el tensor energía-momentum Q üb en la forma canónica 

Q ab = eVaVt + pWaWt, (4.4) 

donde eyp son, respectivamente, la densidad de energía y la presión azimutal del disco. En términos 
de estas cantidades, las condiciones de frontera pueden ser escritas como 

4e 2( *- A) [1 - r* r ] <D, ; = e, (4.5a) 

4e 2(í '- A Vo r <D, z = Pt (4.5b) 
donde, como arriba, todas las cantidades son evaluadas en z = 0 + . 
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Como se puede ver de las expresiones anteriores, el tensor energía-momentum más general com- 
patible con el elemento de linea (4.1), corresponde a una fuente discoidal con densidad de energía y 
presión azimutal diferentes de cero. En acuerdo con esto, en lugar de suponer prescripciones espe- 
cíficas para la densidad de energía y la presión azimutal, se solucionarán las ecuaciones de Einstein 
requiriendo sólo que estas dos cantidades serán diferente de cero en la superficie de un disco con un 
borde interno central. Luego, dada una solución, esta puede usarse para obtener, a partir de las condi- 
ciones de frontera, las correspondientes expresiones para la densidad de energía y la presión azimutal. 
Por consecuencia, las soluciones corresponderán a la clase más general de discos estáticos delgados 
con borde interno central que puedan ser obtenidos solucionando exactamente las ecuaciones de Eins- 
tein. 

De esta manera, debido a que se pretende obtener una solución interpretable como un disco del- 
gado localizado en la hipersuperficie z = 0, con un borde interno central de radio a, sólo es necesario 
imponer que 

í 0 ; 0Sr<«, 
<D,-(r,0 + ) = | (4.6) 
( f(r) ; r>a, 

con f(r) cualquier función arbitraria. Entonces, sólo después de hallar la solución más general, se 
imponen condiciones adicionales dado que se quieren comportamientos físicamente razonables. Así, 
con el fin de obtener un espaciotiempo asintóticamente plano, se requiere que 

lím <£(/-, z) = 0, (4.7a) 

R—too 

límA(r,z) = 0, (4.7b) 

donde R 2 = r 1 + z 2 . También, con el fin de tener regularidad en el eje de simetría, se requiere que 

O(0, z) < oo, (4.8a) 



A(0,z) < oo. (4.8b) 

También se requiere que 

f{r) > 0, (4.9a) 

0 < r<D,,. < 1, (4.9b) 

con tal que la densidad de energía y presión azimutal sean positiva en todas partes. Finalmente, se 
necesita analizar la finitud de la masa total de los discos. Así, con esta intención, primero se toma la 
densidad de masa /j de los discos definida por 

\ = ÍQub - \g<*Q)V a v\ (4.10) 
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donde Q = g" b Q„b y, como antes, todas las cantidades son evaluadas en z = 0 + . Así, usando (4.2) y 
(4.4), se tiene que la densidad de masa se reduce a 

H = e + p. (4.11) 

De otra parte, la masa total de los discos está dada por 

M= f ' iidL= f f fie A - m rdrd<p, (4.12) 

J Jo Ja 

donde dL = e K ~ m rdrdíp es el elemento de área sobre la superficie del disco. Entonces, usando (4.11), 
se puede escribir la masa total como 

M = 2n ^ (e + p)e K ~ m rdr, (4.13) 

donde ya se ha hecho la integración sobre tp. 



4.2. Discos Delgados con Borde Interno Central 

Con el propósito de solucionar las ecuaciones de Einstein en el vacío, en primera instancia se 
tiene que solucionar el problema de valor en la frontera para O. Ahora bien, debido a la paridad de 
las funciones métricas respecto del la superficie del disco, es conveniente escoger un nuevo sistema 
de coordenadas que se adapten naturalmente a la geometría de la fuente deseada. Por tal razón, se 
introducen las coordenadas esferoidales oblatas definidas a través de la relaciones 

r 2 = a 2 (\ +x 2 )(l -y 2 ), (4.14a) 



z = axy, (4.14b) 

donde -oo <jr<ooyO<¡/< 1. Así, cuando x = 0 se tiene que z = 0y0<r<a, mientras que 
cuando y = 0 se tiene que z = 0 y r > a. Además, como x cambia de signo al atravesar la superficie 
y = 0, pero no cambia en valor absoluto, esta coordenada presenta una discontinuidad finita en y = 0. 
Entonces, una función par de x es continua en todas partes pero tiene una derivada normal discontinua 
en y = 0. De otra parte, y es continua en todas partes. En consecuencia, la superficie y = 0 describe 
un disco con un borde interno de radio a, mientras que la superficie x = 0 describe un "hueco vacío" 
con borde en el disco. 

En las coordenadas esferoidales oblatas, el elemento de linea de Weyl (4. 1 ) puede ser escrito como 
ds 2 = - e^dt 2 + a 2 (l+ x 2 )(l - y 2 )e- m d V 2 + «V + y 2 )e 2iA -° } 
de modo que las ecuaciones de Einstein en el vacío se reducen a 

[(1 + x 2 )®, x ], x + [(1 - ¡hfij* = 0- (4-16) 



dx ¿ 



dy~ 



l+x 2 



i-y 2 



(4.15) 
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la ecuación de Laplace en coordenadas esferoidales oblatas, y el sistema sobredeterminado 

A,, = (1- y 2 ) [^(1+^)^-^(1-^)0^ 

- 2y(l + x 2 )<D, x <D.;| l(x 2 + y 2 ), (4.17) 

A, 9 = (l+x^lyd+x^-yd-y 2 )^ 

+ 2x(.l - y 2 )®,®,,] K* 2 + y 2 ), (4.18) 

cuya integrabilidad está garantizada nuevamente por la ecuación (4. 16). 
De otra parte, usando (4.14a) y (4.14b), es fcil ver que 

í <t> -;t (0, y)/ay ; O < r < a, 
O ; (í-,0) = j (4.19) 
I ®, 9 (x, 0)1 ax ; r > a. 

Consecuentemente, como la simetría de reflexión de las soluciones implica que 

0>(-x,y) = <D(x,y), (4.20a) 

<t>,A-x,y) = -®,Ax,y), (4.20b) 
las condiciones (4.6) son equivalentes a 

<5. v (0,¡/) = 0, (4.21a) 

<D, 9 (x,0) = F{x)\ x>0, (4.21b) 

con F(x) una función arbitraria. La solución más general de la ecuación (4.16) con estas condiciones 
de frontera está dada por [2] 

(D(x, y) = J][A 2 „P 2n (y) + B ln Q 2 „(y)]P2„(x), (4.22) 

donde A?„ y B 2 „ son constantes, P 2 „(y) y Q 2 „(y) son los polinomios de Legendre y las funciones de 
Legendre de segunda clase, respectivamente, y p 2n (x) = i~ 2n P 2 „(ix). En consecuencia, todas las solu- 
ciones de las ecuaciones de Einstein en el vacío para espaciotiempos estáticos con cualquier fuente 
axialmente simétrica como la considerada aquí son obtenidas tomando para la función métrica 0(jr, y) 
cualquier elección particular de la solución general anterior, o expresiones obtenidas a partir de estas 
soluciones por medio de operaciones lineales. Ahora, en términos de las coordenadas esferoidales 
oblatas, la condición (4.7a) se escribe como 

lím <D(;t, y) = 0, (4.23) 

X— »oo 

mientras que la condición (4.8a) se escribe como 

<D(x, 1) < co. (4.24) 
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Así, debido al comportamiento de las funciones de Legendre, es claro que no es posible satisfacer las 
condiciones física (4.7) y (4.8) con cualquier escogencia particular de la solución general (4.22). Sin 
embargo, considerando solo el primer término de la serie, 



<b 0 (x,y) = A 0 + B 0 Q 0 (y), 



(4.25) 



se obtiene una solución que es regular para todo y + 1. Entonces, si se toma Aq = 0, esta solución 
puede escribirse como 



■ln 



í + y 



(4.26) 



donde a es una constante arbitraria, de tal modo que una integración directa de (4.17)-(4.18) da 

i-r 



Ao(x, y) : 



■ln 



(4.27) 



donde se ha tomado la constante de integración como cero. Como se puede ver, esta solución es no 
asintóticamente plana ni regular en el eje de simetría. 

Sin embargo, usando (4.5), se obtiene para la densidad de energía y la presión azimutal las expre- 
siones 



(4a/a)x 2<1 



(4.28) 



0. 



(4.29) 



de tal modo que, si a > 0, los discos satisfacen todas las condiciones de energía [39]. Así, para 
cualquier valor de a 2 + |, la densidad de energía crece ilimitadamente, en el infinito o en el borde 
interno del disco; mientras que para alfa a 2 = | la densidad de energía es siempre constante, tal que, 
en cualquier caso, la masa total del disco será infinita. 

De otra parte, aunque las soluciones previas no tienen un comportamiento físico aceptable, es 
posible usar esta como punto de partida para generar una solución nueva y bien comportada. Para este 
fin, se consideran las coordenadas esferoidales oblatas no sólo como funciones de las coordenadas 
cilindricas, (r, z), sino también como paramétricamente dependiente del radio a, 



x = x(r, z\ a), 



(4.30a) 



y = y(r,z;á). 

Asimismo, suponiendo también que í> depende paramétricamente de a, 

O = 0(r, z; a), 

se puede obtener una nueva familia de soluciones aplicando la operación lineal 



í>„ +1 (r,z;a) ; 



da 



(4.30b) 



(4.31) 



(4.32) 
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donde n es un entero, n > 0. 

Entonces, a partir de la "solución semilla" OoC*, í/)> es posible generar, por medio del procedi- 
miento previo, una nueva familia de soluciones que pueden ser escritas en la forma simple 

para n > 1, donde los FJx, y) son funciones polinómicas cuyo mayor grado es 4n - 4, de las cuales 
los tres primeros son, 



1, 

x 4 + 3x 2 (l - y 2 ) - y 2 , 



F 3 = 3 x 6 (3 - 5?/ 2 ) + 5x 4 (6í/ 4 - Uy 2 + 3) 
- x 2 y 2 (3y 4 - 3ly 2 + 30) - y\y 2 - 3), 



pero todos ellos puede obtenerse fácilmente por medio de (4.32). Así, es fácil ver que 

lím<D„(x,¡/) = 0, (4.34) 

x— >eo 

y que 

<D„(x, 1) < oo, (4.35) 

están en completo acuerdo con las condiciones (4.7a) y (4.8a). 

Ahora, con el fin de obtener las correspondientes funciones métricas A„(r, z; a), se hace la integra- 



A„(r, z; a) = A„(x, V) = J A .y( x > y) d V> ( 4 - 36 ) 

tomando A„(x, 1) = 0 para garantizar la regularidad en ele eje. Así, poniendo (4.33) en (4.18), las 
soluciones pueden ser escritas en la forma simple 

. . , a 2 (2n-2)\(y 2 - l)A„(x,y) 

K(X > y) = A^ { x 2 + y 2 f- ' (437) 

para n > 1, donde los A„(x, y) son funciones polinómicas, cuyo mayor grado es 8n - 2, de los cuales 
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se presentan aquí sólo los tres primeros 

A, = x\9y 2 - 1) + 2x 2 y 2 (y 2 + 3) + y\y 2 - 1), 

A 2 = 2x l2 (9y 2 - l)-4x 10 (51¡/*-41¡/ 2 + 2) 

+ x g (735if - \2A\y 4 + A\9y 2 - 9) - A 2 (132¡/ 6 

- 1644¡/* + 1604y 2 - 252) + x 4 y 4 (84y 6 - 384¡/ 4 
+ 1266y 2 - 630) + 4x 2 if(6if + 6y 4 - 39y 2 + 63) 

+ 3í/V + y 4 + y 2 - 3), 

A 3 = 3x 16 (1225í/ 6 - 1275¡/ 4 + 315¡/ 2 -9) 

- 24x 14 (980i/ 8 - 2095¡/ 6 + 1205i/ 4 - 189¡/ 2 + 3) 
+ 2x 12 (24255i/° - 89475y 8 + 98472¡/ 6 - 363 16¡/ 4 
+ 3473y 2 - 25) - 12x 10 ¡/ 2 (1835í/ 10 - 16665í/ 8 
+ 34716y 6 - 25292í/ 4 + 6001i/ 2 - 275) 
+ 6x s y\900y 10 - 1 1946¡/ 8 + 50563í/ 6 - 69397í/ 4 
+ 33365¡r - 4125) + 8xV(125¡/ 10 + 926¡/ 8 

- 9079y 6 + 24639¡/ 4 - 22290i/ 2 + 5775) 
+ 6*y (55í/ 10 + 29¡/ 8 + 764¡/ 6 - 4808¡/ 4 + 8469t/ 2 

- 4125) + n^y'^Sí/ 10 + 5if + SOy 4 - 30ly 2 
+ 275) + y 12 (5y 10 + 5if + 5y 6 + y 4 + 34y 2 - 50), 

todos ellos pueden obtenerse como un resultado del cálculo de la integral (4.36) 
tiene que 

lím A„(x, y) = 0, 

y que 

A„(x, 1) = 0, 

en completo acuerdo con las condiciones (4.7b) y (4.8b). 

4.2.1. Comportamiento de las Soluciones 

Como se puede ver de las expresiones en la sección previa, usando (4.32) se generó una familia 
infinita de soluciones asintóticamente planas, todas con un comportamiento regular en el eje de sime- 
tría. De otra parte, parece ser que existe una singularidad en el borde interno de los discos, cuando 
x = 0 and y = 0. Sin embargo, es fácil ver que esta aparente singularidad es sólo una singularidad 
coordenada. Ahora bien , como es ampliamente conocido de [84, 95], la caracterización invariante de 
la curvatura puede darse en términos de 14 escalares construidos a partir del tensor de Riemann, S„t„/, 
y el tensor métrico, g a ¡,. Así que, en orden a determinar la naturaleza de las aparentes singularidades, 
se tiene que calcular todos los invariantes de curvatura para la familia de soluciones. 



. Acordemente, se 
(4.38) 

(4.39) 
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Ahora, para cualquier solución de las ecuaciones de Einstein, los invariantes de curvatura son los 
10 componentes nulos del tensor de Ricci, R ü ¡, = 0, más los cuatro escalares 

ti/" nclbcd n 

A/ - K K abal , 

ti/" nub nklcdn 

A// - K k ,K K abííh 



III 



donde g = det g„¡, y é' bcd es el símbolo de Levi-Civita. Además, para cualquier solución de Weyl los 
dos últimos invariantes se anulan idénticamente, de modo que sólo es necesario calcular TC; y "Ka. 

Usando las expresiones (4.33) y (4.37) para 0„(x, y) y A„(x, y), se pueden tener los escalares de 
curvatura como 

16a 2 e^-^N In (x,y) 
% " = ~ a^(x 2 + y 2 y 2 » ' (44 ° a) 

48a^-^N, ln (x,y) 
%1 " ~ a^(x 2 + y 2 )^ ' (440b) 

donde N¡„(x,y) y N¡i„(x,y) son funciones polinómicas, cuyo máximo orden es 24íi - 6 y 32íi - 9, 
respectivamente, los cuales no se escriben debido a su gran extensión, pero que todos se anulan en el 
borde interno de los discos, 

A' / „(0,0) = A'„„(0,0) = 0. (4.41) 

Además, es fácil ver que, en cualquier vecindad alrededor de (0,0), la diferencia entre <¡>„(x, y) y 
A„(x, y) se comporta como 

a 2 

1>„ - A„ » - — — , (4.42) 

a 2 "(x 2 + y 2 ) 2 " 



de tal modo que 



lím — — = 0, (4.43a) 

(x,y)^(0,0) (x 2 + y 2 ) [2 " 

g6(Q„-A„) 

lím t--, — = °' ( 4 - 43b ) 

U»)-KP,0) (x 2 + y 2 ) 16 " 

y los límites existen, cualquiera que sea la trayectoria escogida para aproximarse al punto (0, 0). 
Acordemente, se tiene que 

lím %„(x,y) = 0, (4.44a) 

(*,jí)->(0,0) 

lím %u„(x,y) = 0, (4.44b) 
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y entonces la curvatura es regular en el borde interno de los discos. 

Ahora, con el fin de analizar el comportamiento físico de las fuentes, se calculará la densidad de 
energía y la presión azimutal para estas familias de discos. Así, usando (4.14a) y (4.14b), se puede 
ver que 

Vi + x 2 

®Jx,0), r> a, (4.45) 



O r (r,0) : 

y, usando (4.33), es fácil probar que 

<D„v(x,0) = 0, n>\. (4.46) 

Entonces, de (4.5b), se puede ver que 

Pn = 0. (4.47) 

Esto es, todos los discos de la familia tienen presión azimutal cero. 

De otro lado, usando las ecuaciones (4.5a), (4.19) y (4.33), la densidad superficial de energía de 
los discos puede ser escrita como 

AaE„(x) _____ f a\2n-2)\B„{x) 
a" +l x- 



e » W = ^„ +lv 2„ +1 ex P i 2 2 "- l a 2 "x 4 " í ' (4 ' 48) 



donde x > 0 y las cantidades E„(x) positivo definidos de grado 2k, con k = (n - l)/2 para n impar y 
k = n/2 para n para, de los cuales a continuación se escriben sólo los tres primeros términos, 

Ei{x) = 1, 
E 2 (x) = x 2 + 3, 
E 3 (x) = 3(x 2 + 5), 

mientras que los B„{x) son polinomios positivo definidos de grado Ak, con k = (n - l)/2 para n impar 
y k = n/2 para n impar, de cuyos tres primeros tres términos están dados por 

Bi(x) = 1, 

B 2 (x) = 2x 4 + Sx 2 + 9, 

B 3 (x) = 27x 4 + 72x 2 + 50. 

De las expresiones anteriores se puede ver que, tomando a > 0, la densidad de energía de los 
discos será positiva siempre, 

eM > 0. (4.49) 

Así, como la presión azimutal es cero, se tiene una familia infinita de discos de polvo que están en 
completo acuerdo con todas las condiciones de energía. También es fácil que, para cualquier valor de 
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n, se tiene que 

e„(0) = 0, (4.50a) 

líme„(x) = 0. (4.50b) 

Esto es, la densidad de energía de los discos es cero en su borde interno y se anula en el infinito. 
Además, como la presión azimutal, la densidad de masa de los discos se reduce a su densidad de 
energía, 

H n (x) = e n (x), (4.51) 

tal que su comportamiento es el mismo que la densidad de energía. 

Ahora, con el fin de mostrar el comportamiento de la densidad de energía, se gráfica la densidad 
superficial de energía adimensional e„ = ae„ como funciones de la coordenada radial adimensional 
f = r/a. Así, en la figura 4.1, se gráfica e„ como una función de f para los primeros tres términos de 
la familia, con n = 1, 2 y 3, para diferentes valores de los parámetros á„ = a/a". Entonces, par cada 
valor de n, se tiene á„ = 0,5, 1, 1,5, 2, 2,5, 3, 3,5 y 4. La primera curva de la izquierda corresponde a 
á„ = 0,5, mientras que la última curva en el lado derecho corresponde a á„ = 4. Tal como se puede 
ver, en todos los casos la densidad de energía superficial es siempre positiva, con un máximo cerca 
al borde interno de los discos, y luego decrece rápidamente como f crece. Se puede ver que, para 
un valor fijo de n, conforme los valores de á„ crecen, los valores de los máximos disminuyen y se 
mueven hacia los valores en que f crece. El mismo comportamiento se observa para un valor fijo de 
á„ incrementando n. 

De otra parte usando, (4.13), (4.14a) y (4.51), la masa de los discos puede expresarse como 



•f 



M„ = 2na \ f„(x)dx, (4.52) 
donde 

/„(*) = xe- K ^ m e n (x). (4.53) 
Así, de (4.37), (5.2) y (4.53), es fácil ver que 

/n+lW £„+¡(x) . 

hm = hm , (4.54a) 

fn(x) "» e n (x) 

x-k*> ax ¿ h„(x) 

[ \ ; n = 2k+\, 



(4.54c) 



0 ; n = 2k + 2, 



con k > 0. Por consecuencia, mediante el criterio de comparación de los límites para integrales 
impropias, la convergencia de M n+ ¡ está garantizada si M„ es convergente, y entonces sólo es necesario 
probar la convergencia de M¡ . En efecto, mediante un simple cálculo se tiene 

M¡ = 2n V2aa T( 1 /4), (4.55) 
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n= 1 




1 15 2 2.5 3 



Figura 4.1: Densidad superficial de energía e„ en función de f para los tres primeros discos de la 
familia, con á„ = 0,5, 1, 1,5, 2, 2,5, 3, 3,5 y 4. Para cada valor de n, la primera curva a la izquierda 
corresponde aí = 0,5, mientras que la última curva a la derecha corresponde aá = 4. 
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lo que garantiza así la convergencia de las integrales de masa (4.52). Consecuentemente, aunque los 
discos son de extensión infinita, todos tienen masa finita. 

En conclusión, se ha presentado una familia infinita de soluciones de las ecuaciones de Einstein 
en el vacío asintóticamente planas y regulares en todas partes. Estas soluciones describen una familia 
infinita de discos delgados con un borde interno central constituidos por materia de polvo, cuya densi- 
dad de energía es siempre positiva y bien comportada, de tal modo que su tensor energía-momentum 
está en completo acuerdo con todas las condiciones de energía. Además, aunque los discos son de 
extensión infinita, todos tienen masa finita. 

Ahora, como las funciones métricas de las soluciones son calculadas explícitamente, estas cons- 
tituyen la primera familia de soluciones exactas explícitas completamente integradas para tal clase 
de fuentes de discos delgados. Además, su relativa simplicidad cuando son expresadas en coordena- 
das esferoidales oblatas, permite estudiar muy fácilmente sus diferentes aspectos dinámicos, como el 
movimiento de partículas dentro y fuera de los discos y la estabilidad de las órbitas. 

Aparte de su importancia como una nueva familia de soluciones exactas y soluciones explícitas de 
las ecuaciones de Einstein en el vacío, la principal importancia de esta familia de soluciones es que ell 
puede ser fácilmente superpuesta con la solución de Schwarzschild con el fin de obtener la descripción 
de un sistema binario compuesto por un disco delgado alrededor de un agujero negro central. En 
efecto, la superposición del primer miembro de esta familia con un agujero negro de Swarzschild 
fue obtenida y presentada en in [31]; en la próxima sección se presentará un análisis detallado de 
la superposición de la familia completa. Consecuentemente, como las familias presentadas aquí, su 
superposición con el agujero negro de Swarszchild será la primera familia de soluciones exactas 
integradas explícitamente para esta superposición de fuentes. 

Ahora, como en la teoría newtoniana el potencial gravitacional está dado por la solución del pro- 
blema de valores en la frontera para la ecuación de Laplace, se puede considerar Q„(r, z\ a) como una 
familia de potenciales gravitacionales newtonianos de fuentes discoidales con borde interno, cuyas 
densidades de masa newtoniana están dadas por 

2aE„(x) 

= a ,+i x 2„ + i ' ( 4 - 56 ) 

claramente divergente en el borde de los discos. Por consecuencia, se puede concluir que no existen 
soluciones en la teoría newtoniana que describan el campo gravitacional de un disco delgado con 
borde interno, mientras que esta clase de fuentes puede ser propiamente descrita por medio soluciones 
regulares y asintóticamente planas, en el marco de la teoría de la gravitación ofrecida por la relatividad 
general. 



4.3. Discos de Polvo Contrarrotantes Alrededor de un Agujero 
Negro de Schwarzschild 

A continuación se presenta la primera familia de soluciones exactas, completamente integrada 
explícitamente, de las ecuaciones de Einstein correspondiente a la superposición de un disco de polvo 
contrarrotante con un agujero negro. La solución obtenida representa un disco anular delgado infinito 
(un disco con un borde interno) alrededor de un agujero negro de Schwarzschild. La masa del disco es 
finita y el tensor de energía-momentum cumple con todas las condiciones de energía. Además, la masa 
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total del disco en presencia del agujero negro es menor que cuando el disco está solo. La solución 
puede interpretarse también como si describiera un disco delgado constituido por dos fluidos de polvo 
contrarrotantes alrededor de un agujero negro que satisfacen también las condiciones de energía. 

Como es sabido, las ecuaciones de campo de Einstein para la métrica de Weyl-Lewys-Papapetrou 
(2.5) en las coordinates x" = (f, tp, r, z), son 

O, T + -<D,. + <í>- = 0, (4.57a) 
r 

A ,. = r (Of r - <D 2 ), A, ; = 2rO, r O, ; . (4.57b) 

Ahora se procederá a superponer los discos de polvo con borde interno central <é„, estudiados en 
secciones precedentes, alrededor de un agujero negro ip. Ahora bien, la ecuación (4.57a) es sólo la 
ecuación de Laplace en coordenaas cilindricas. La ecuación de Laplace es lineal, de modo que si <p n y 
(A son soluciones, lo será también la superposición <E>„ = <f>„ + \¡i . 

No obstante, se puede ver de (4.57b) que, la función métrica A[0„] no puede ser superpuesta dado 
que esta no es lineal. Sin embargo, se puede demostrar que esta obedece la relación 

A[<f>„ + #] = A[0„] + A[<A] + 2A[<A„,(A], (4.58) 

donde 



A[<A„, l/f] = J T [(|/V0„.,- " ^,z<f>,;z) dr + (>l>,r<t>n¿ + <l< ' An,r)d¿[, 



(4.59) 



y el primer modelo de discos de polvo con borde interno central es 



Mx,y,a)= " y (4.60) 
a(x 2 + y 2 ) 

siendo o- una constante. Ahora, usando (4.58) y (4.5) se puede calcular el tensor energía-momentum 
del sistema combinado <5„ = (f>„ + tp. Así, la densidad superficial de energía, la presión azimutal y la 
"densidad newtoniana efectiva" del sistema compuesto es: 

e = 4e 2 *"- 2A[ '*" +w (l-^„ + i/ í ) J .}(0„ + iA)_ ;: , (4.61a) 
p = Are 1 "'- 2 ^^ + f\ r (<t> n + <A), ; , (4.61b) 
p = e + p = 4e 2 *"- 2A[ *" +,íl (0„ + iA), ; , (4.61c) 

donde A[(f>„ + ijj] es dado por (4.58) y todas las cantidades son evaluadas en z = 0 + . 

De otra parte, la función métrica \p correspondiente a un agujero negro de Schwarzschild no ro- 
tante viene dado por [84] 



,= Il„{|^f, (4.62) 

donde ( y r¡ son las coordenadas esferoidales prolatas, relacionadas con las coordenadas cilindricas 
mediante 

r 1 = m 2 (f - 1)(1 - í? 2 ), z = mfr, (4.63) 
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donde 1 < f < oo y — 1 < 37 < 1, con m una constante positiva. El otro potencial métrico es 

hW = \\n{p^\. (4.64) 



2 (C 

Resulta de gran utilidad reexpresar las cantidades (4.62) y (4.64) en una forma práctica, que como se 
verá pronto, facilitará los cálculos en el proceso de superposición. Con esto en mente, se definen las 
cantidades 

/íj = -m-z + Ri (4.65a) 
p 2 = m - z + R 2t (4.65b) 



R 2 ¡ = r 2 + (m + zf (4.66a) 
R\ = r 2 + (m-z) 2 . (4.66b) 

A partir de aquí, se obtienen las siguientes identidades útiles: 

Rt + R*> Rt — R*> 

2m 2m 

entonces, con (4.65b) y (4.67) se obtiene de (4.62) y (4.64) respectivamente 

* = 5 ln {£}' (468a) 

AfiH = f + £*? J . («Sb) 

que serán las fórmulas a utilizar en la supeiposisción de los discos alrededor del agujero negro. 

Ahora, con el fin de obtener el término mixto A[<p n , 41] en (4.58) en la superposición, se introduce 
la relación entre las coordenadas esferoidales prolatas y las cilindricas mediante 

R + R' R-R* 
2a 2a¡ 



R = Vr 2 + (z + ai) 2 , R* = ^r 2 + (z - ai) 2 . (4.70) 

Luego, se definen las cantidades auxiliares 

fi = R-(z + aí), m* =R* ~(z-ai). (4.71) 

Así, introduciendo las anteriores expresiones en (4.60), se puede expresar la función (f>¡, a través 
de la cantidad /i, como 

0i=-~ (lníí + ln/O. (4.72) 
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Entonces, usando (4.68a) y (4.72), se obtiene para A[</>i,<//]: 

a d 



A[<A,,^] = A 



2 ¡J. 2 



2 da 



(InjU + ln//) 



ahora, unos pocos pasos de álgebra elemental dan 



a 8 
~Adli 



ln 



a(l-y){A l -A 2 ) 
2(x 2 + y 2 ) ' 



donde 



¿i 

¿2 



ay(l - y)(x 2 + 1) - mx(\ + y){( - 1)(1 - rj) 

[ax + m(l - £ - r¡)] 2 + a 2 (l - y) 2 
ay(l - y)(x 2 + 1) - mx(l + y){( + 1)(1 - r¡) 



[ax - m(l + <f - ij)] 2 + a\\ - y) 2 



y se ha hecho uso de las propiedades 



F[a(f>] 
F[a(f>i;(f>2] 



d_ 
da 



F[Ma);<f>2] 



a 2 F[cf>], 

aF[(p\;(j>2\, 

F[^;0 3 ]+F[^ 2 ;^], 

F[— ^;0 2 ] 



(4.73) 



(4.74) 



(4.75a) 
(4.75b) 
(4.75c) 

(4.75d) 



con a un parámetro arbitrario y a una constante. 

Así, si se toma fa, el primer miembro de la familia de discos con borde interno central, se puede 
superponer este a un agujero negro de Schwarzschild usando (4.74) en (4.58). 

De otra parte, se sabe de (4.32), que si es posible obtener el n + 1-miembro de la familia de discos, 
0„+i, a través de la fórmula de recurrencia dada por 



¡p„ +l (x,y;a) ■■ 



da 



4>„(x,y;a), 



entonces, es posible construir una fórmula de recurrencia para la superposición . 
O n+ i = 0„ + i + ijj. Consecuentemente, para el otro potencial métrico se tiene que 



A[0„ +1 + i/i] 



A[0„ +1 ]+A[<A] + 2A[<fc +1 ,iA] 

d 3 
A[— « + A[<A] + 2A[— <f>„M 



da 



da 



da da 



(4.76) 
+ ijj, dado que 



(4.77) 



donde nuevamente se ha hecho uso de las propiedades (4.75). Así entonces, a partir de la super- 
posición del primer miembro de la familia de funciones métricas h[t¡>\ + ¡A], se puede obtener la 
superposición correspondiente al íi-miembro A[^„ + (A] a través de derivaciones sucesivas. 
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De otro lado, se puede probar que la solución para el término mixto de la superposición, A[(f>„, ifr], 
evaluado en y = r¡ = 0 es A[(p n , i//\ y=n =o = 0. De igual modo, las derivadas de las funciones métricas 
<Pn, r y (A,z se anulan en y = n = 0. Por consecuencia, la densidad superficial de energía, la presión 
azimutal y la densidad newtoniana efectiva del sistema superpuesto pueden ser escritos como 

e„ = e 2 ^ 2A ^(l-r^^„, (4.78a) 
p„ = nA, r e 2, *- 2AW e„, (4.78b) 
p„ = e 2 *- 2A M £ „, (4.78c) 



s n = §^„, (4.79a) 



(<r+D 2 
i! 

(<r+D 2 



P- = 77^ e n, (4.79b) 



con p„ = e„ + p„, siendo 



AaE n (x) <_a 2 (2n-2y.B„(x)) 
eÁX) ~ a" + 'x 2 " +l BXP \ l^a 2 "^» ) ' 



donde x > 0 y los E„(x) son polinomios positivo definidos de grado 2k, con k = (n - l)/2 para n 
impar, y k = n/2 para n par; de los cuales los tres primeros tres son expuestos a continuación: 

E¿x) = 1, 
E 2 (x) = x 2 + 3, 
E 3 (x) = 3(3x 2 + 5), 

mientras que los B„(x) son polinomios positivo definidos de grado 4k, con k = (n- l)/2 para n impar, 
y k = n/2 para n par; los primeros tres de ellos son dados por 

Bito = 1, 

B 2 (x) = 2x 4 + 8x 2 + 9, 

B 3 (x) = 27x 4 + 72x 2 + 50. 

Ahora, como e>0yf>l, eyp serán siempre positivos. Consecuentemente, el tensor energía- 
momentum estará en completo acuerdo con las condiciones de energía fuerte y débil. De otra parte, 
con el fin de satisfacer la condición de energía dominante se requiere que p < e, lo cual implica que 
C > 2 y a > V3m. Además, se tiene que p < e, y entonces la masa total del disco en presencia del 
agujero negro es menor que la masa del disco solo. 
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El tensor energía-momentum puede ser interpretado como la superposición de dos fluidos contra- 
rrotantes. Con el fin de estudiar esto, se escribe Q" 1 ' como [22] 

Q" b = e + U'l U b + + e_ U'L Ut, (4. 80) 

donde 



% (4.81) 



son las densidades energía de lso dos fluidos contrarrotantes. Las velocidades de contarrotación son 
dadas por [22] 

V" ± UW" 

U: = . (4.82) 

donde 

U 1 = - = < 1 (4.83) 
e £ - 1 

es la velocdad tangencial de contrarrotación. Así, se tienen dos fluidos de polvo con densidades super- 
ficiales de energías iguales. Ahora, como e ± > 0, los dos discos de polvo están en completo acuerdo 
con todas las condiciones de energía. Como se puede ver, las soluciones anteriores presentan algu- 
nas propiedades interesante. Primero, las fuentes de materia asociadas presentan un comportamiento 
muy razonable y su tensor energía-momentum está en completo acuerdo con las condiciones de ener- 
gía. Además, su relativa simplicidad cuando son expresadas en términos de coordenadas esferoidales 
prolatas y oblatas, permite el fácil estudio de aspectos dinámicos, como por ejemplo el movimien- 
to de las partículas dentro y fuera del disco, la estabilidad de las órbitas y la posible existencia de 
singularidades. 

Ahora, con el fin de mostrar un análisis del comportamiento de las densidades superficial de 
energía de la superposición, se graficará la densidad superficial de energía adimensional e„ = as„ 
en función de la coordenada radial adimensional f = r/a. Así, en la Figura 4.2, Figura 4.3 y en la 
Figura 4.4 se gráfica e„ en función de f para los tres primeros miembros de la familia de soluciones 
correspondientes a la superposición, con n = 1, 2 y 3 respectivamente, para diferentes valores de 
los parámetros m y a„ = a/a". Así, para cada valor de n, se toma m = V3a/3, V3~a/6, y V3a/9 y 
a n = 1 ,2,. . ., y 9. La primera gráfica de arriba hacia abajo en Figura 4.2 corresponde a m = V3a/3. 
La primera de las curvas, de arriba hacia abajo, en esta gráfica corresponde ¿1 = 1, mientras que la 
última curva, corresponde a ág = 9. Como se puede ver, en todos los casos la densidad superficial de 
energía es siempre positiva, con un máximo cerca del borde de los discos, y luego decrece rápidamente 
conforme f decrece. Se puede ver que, para un valor fijo de n y m, conforme los valores de á„ crecen, 
el valor del máximo disminuye y se desplaza hacia donde crecen los valores de f. De otra parte, para 
un valor fijo de á„ y m el máximo crece y se desplaza en dirección de crecimiento de f, conforme n 
crece. Mientras que para n y á„ fijos, los máximos crecen como m decrece. 

Posteriormente se gráfica la presión azimutal p„ = ap„ como función de la coordenada radial 
adimensional f = r/a. Así, en la Figura 4.5, Figura 4.6 y Figura 4.7 se gráfica p„ como una función de 
f para los tres primeros miembros de la familia de la superposición, con n = 1, 2 y 3 respectivamente, 
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para diferentes valores de los parámetros m y á„ = a/a". Así, para cada valor de n, se toma m = 
V3a/3, VIa/6, y V3a/9 y á„ = 1 , 2,. . . , y 9. La primera gráfica de arriba hacia abajo en la Figura 
4.5 corresponde a m = V3a/3. La primera curva de arriba hacia abajo en esta gráfica corresponde a 
¿1 = 1, mientras que la última en esta misma gráfica corresponde a = 9. Como se puede ver, en 
todos los casos la presión azimutal es siempre positiva, con un máximo cerca del borde interno de 
los discos, y luego decrece rápidamente conforme f crece. Se puede notar también que, para valores 
fijos de n y m, conforme los valores de á„ crecen, los valores máximos de la presión disminuyen y se 
desplazan en la dirección de crecimiento de f. De otra parte, para valores fijos de á„ y m el máximo 
de la presión decrece y se desplaza en dirección de crecimiento de f conforme n crece. Mientras que 
para n y á„ y, el máximo decrece conforme m decrece. 

Finalmente se gráfica la densidad newtoniana efectiva adimensional p„ = ap„ en función de la 
coordenada radial adimensional f = r/a. Así, en la Figura 4.8, Figura 4.9 y Figura 4.10 se gráfica 
p„ en función de f para los tres primeros miembros de la familia de soluciones correspondientes a 
la superposición, con n = 1, 2 y 3 respectivamente, para diferentes valores de los parámetros m y 
a„ = a/a". Así, para cada valor de n, se toma m = V3a/3, V3a/6, y V3fl/9 y a„ = 1 , 2,. . . , y 9. 
La primera gráfica de arriba hacia abajo en Figura 4.8 corresponde a m = V3a/3. La primera curva 
en esta gráfica corresponde a á¡ = 1, mientras que la última corresponde a ág = 9. Como se puede 
ver, en todos los casos la densidad newtoniana efectiva es siempre positiva, con un máximo cerca del 
borde interno de los discos y decreciendo rápidamente conforme f crece. También se puede ver que, 
para un valor fijo de n y m, conforme los valores de á„ crecen, el valor del máximo de la densidad 
newtoniana disminuyen y se desplazan en la dirección de crecimiento de f. Mientras que para «ya, 
fijos, el máximo crece como m decrece. 

De modo que, se ha presentado un formalismo para superponer discos de polvo con un borde 
interno central alrededor de un agujero negro de Schwarzschild. la superposición se obtuvo usando 
la solución de las ecuaciones de campo de Einstein correspondientes a discos de polvo con un borde 
interno central del trabajo [25]. El formalismo para la superposición fue inspirado en el trabajo de 
Letelier and Oliveira [56]. 

Las densidades de energía del sistema superpuesto son siempre funciones positivas de la coorde- 
nada radial, igual a cero en el borde interno central del disco, con máximo en un valor finito de los 
radios y cayendo a cero en el infinito. Los discos tienen masa finita y su tensor energía-momentum 
está de acuerdo con todas las condiciones de energía. Ahora bien, la presión azimutal y la densidad 
newtoniana efectiva de los discos tienen comportamiento similar al de la densidad de energía, así que 
el análisis previo es aplicable también a estas cantidades. De otra parte, la velocidad de contrarrota- 
ción de las partículas es siempre positiva y menor que la velocidad de la luz, cayendo a cero en el 
infinito. 
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Figura 4.2: Densidad superficial de energía s = as en función de f = r/a para el primer modelo de 
discos, n = 1, con borde interno central de radio a alrededor de un agujero negro de masa m = 
primera figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica 
s{f) para 1 < r < 3 con diferentes valores del parámetro á„ = a/a". Primero se toma a¡ = 1, la 
primera curva para cada valor de m, luego á 2 = 2,. . . ,y á 9 = 9. 
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Figura 4.3: Densidad superficial de energía s = as en función de f = r/a para el segundo modelo de 
discos, n = 2, con borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, 
primera figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica 
s(f) para 1 < f < 3 con diferentes valores del parámetro á„ = a/a". Primero se toma a¡ = 1, la 
primera curva para cada valor de m, luego á 2 = 2,. . . ,y á 9 = 9. 
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Figura 4.4: Densidad superficial de energía s = as en función de f = r/a para el tercer modelo de 
discos, n = 3, con borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, 
primera figura de arriba hacia abajo, luego m = V3«/6 y finalmente V3a/9. En cada caso se gráfica 
s(r) para 1 < r < 3 con diferentes valores del parámetro á„ = a/a". Primero se toma a¡ = 1, la 
primera curva para cada valor de m, luego á 2 = 2,. . . ,y á 9 = 9. 
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Figura 4.5: Presión azimutal p = ap en función de f = rja para el primer modelo de discos, n = 1, 
con borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, primera 
figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica p{f) para 
1 < f < 3 con diferentes valores del parámetro a„ = a/a". Primero se toma &j = 1, la primera curva 
para cada valor de m, luego á 2 = 2,. . . , y á 9 = 9. 
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Figura 4.6: Presión azimutal p = ap en función de f = r/a para el segundo modelo de discos, n = 1, 
con borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, primera 
figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica p{f) para 
1 < f < 3 con diferentes valores del parámetro a„ = a/a". Primero se toma &j = 1, la primera curva 
para cada valor de m, luego á 2 = 2,. . . , y á 9 = 9. 
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Figura 4.7: Presión azimutal p = ap en función de f = r/a para el tercer modelo de discos, n = 3, con 
borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, primera figura de 
arriba hacia abajo, luego m = V3a/6 y finalmente V3ti/9. En cada caso se gráfica p{f) para 1 < f < 3 
con diferentes valores del parámetro á„ = a/a". Primero se toma ái = 1, la primera curva para cada 
valor de m, luego a 2 = 2,. . . , y a 9 = 9. 
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Figura 4.8: Densidad newtoniana efectiva p = ap en función de f = rja para el primer modelo de 
discos, n = 1, con borde interno central de radio a alrededor de un agujero negro de masa m = 
primera figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica 
p(f) para 1 < f < 3 con diferentes valores del parámetro a n = a /a". Primero se toma á¡ = 1, la 
primera curva para cada valor de m, luego a 2 = 2,. . . , y á 9 = 9. 
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Figura 4.9: Densidad newtoniana efectiva p = ap en función de f = r/a para el segundo modelo de 
discos, n = 2, con borde interno central de radio a alrededor de un agujero negro de masa m = V3a/3, 
primera figura de arriba hacia abajo, luego m = V3a/6 y finalmente V3a/9. En cada caso se gráfica 
p(f) para 1 < r < 3 con diferentes valores del parámetro a n = a /a". Primero se toma á¡ = 1, la 
primera curva para cada valor de m, luego a 2 = 2,. . . , y á 9 = 9. 
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Figura 4.10: Densidad newtoniana efectiva p = ap en función de f = r/a para el tercer modelo de 
discos, n = 3, con borde interno central de radio a alrededor de un agujero negro de masa m = 
primera figura de arriba hacia abajo, luego m = V3«/6 y finalmente V3a/9. En cada caso se gráfica 
p(f) para 1 < f < 3 con diferentes valores del parámetro a n = a /a". Primero se toma á¡ = 1, la 
primera curva para cada valor de m, luego a 2 = 2,. . . , y á 9 = 9. 
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Superposición de Discos de Polvo 
Contrarrotantes Alrededor de un Agujero 
Negro de Schwarzschild Cargada 



Indice 

5.1. Introducción 63 

5.2. Superposición Cargada a través de los Potenciales Complejos de Ernst 63 



5.1. Introducción 

continuación se presenta, lo que constituye hasta donde se sabe, la primera familia de solucio- 
nes exactas completamente integradas explícitamente de las ecuaciones de Einstein-Maxwell 
correspondientes a la superposición de un disco de polvo contrarrotante con un agujero negro. Las 
soluciones obtenidas representan la superposición cargada de un disco infinito delgado anular (un 
disco con un borde interno) alrededor de un agujero negro de Schwarzschild. La masa total del disco 
es finita y el tensor energía-momentum está de acuerdo con todas las condiciones de energía. Además 
los estados asintóticos de la superposición cargada se reducen a la superposición sin carga eléctrica 
descrita en el capítulo anterior. La solución puede interpretarse como la descripción de un disco cons- 
tituido por fluidos provistos de carga eléctrica contrarrotando alrededor de un agujero negro y que, 
están también en completo acuerdo con las condiciones de energía. 



5.2. Superposición Cargada a través de los Potenciales Comple- 
jos de Ernst 

Mediante el uso del "método de los potenciales complejos de Ernst [14, 15] ahora se obtendrá una 
solución de las ecuaciones de Einstein-Maxwell para un espaciotiempo estático axialmente simétrico. 
Como es bien sabido, soluciones de electrovacío de las ecuaciones de Ernst 

(£f-/? 2 )V 2 f = 2TVf-Vf, (5.1) 

conducen a soluciones estacioarias, axialmente simétricas de las ecuaciones de Einstein-Maxwell. La 
constante /3 se define como 0 1 = 1 - qq* con q = p + iq, tal que p y q son reales. Asterisco denota 
complejo conjugado. En el caso particular aquí discutido, un espacio-tiempo estático, ¿¡ = f , las 
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funciones métricas en (2.5) y el potencial electroagnético pueden obtenerse a partir de los potenciales 
complejos de Ernst como sigue: 

f2_ fí 2 

(5.2a) 



(f+D 2 



A '* = & +f)te-P 2 f [fA ^ + 1) ~ ^ X{1 ~ l/) ~ 2y(x2 + "°' '• (5 ' 2b) 



(x + 1)6 

A . 9 = (X 2 ; y2 )( f_ /í 2 ) 2 &y(x 2 + 1) - &d - r) + 2*1 - ,/x..t.\. isa-i 



1 = (íTTyt ' (5 ' 2d) 

* " <jTT)¡ ' (5 - 2e) 

A, = ^, (5.2f) 

donde se ha hecho uso de las coordenadas oblatas (x, y ; a). Naturalmente, si/J = 1 se tienen solucio- 
nes de vacío de las ecuaciones de Ernst. Ahora bien, suponiéndose que se tiene una solución de las 
ecuaciones de vacío de Ernst ¿¡o entonces la transformación 

(5.3) 

con f} 2 < 1 conduce a una solución f de las ecuaciones (5.1). Una simple e interesante escogencia de 
¿¡o en (5.3) es 

f 0 = -coth(£), (5.4) 

donde (f> es una solución de la ecuación de Laplace. Después de la transformación f = -yScoth(i^) los 
potenciales métricos y el potencial eléctrico son 



2/3e* 



«2*(1 -/?)-(!+£)' 



(5.5a) 



A -' = , U 2 V i [x(x 2 + Dft - x(\ - y 2 )^ - 2y(x 2 + Dfr^J, (5.5b) 



{x 2 + y 2 ) 
(x 2 + 1) 
(x 2 +y 2 ) 



[y(x 2 + 1)^ - y(l - y 2 )<j> 2 tJ + 2x(l - y 2 )<t>, x <t>J, (5.5c) 



Vi -Pie 24, - 1) 

A > - Ai-fí-w (5 - 5d) 

Ahora bien, si en (5.4) se usa como solución de las ecuaciones de Einstein en el vacío correspon- 
diente a la superposición de un disco con borde interno central alrededor de un agujero negro de 
Schwarzschild estudiada en el capítulo anterior (í>„, A[0„]) 

<A ^ <D„ = c/>„ + i/f, (5.6) 
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la solución correspondiente a las ecuaciones de Ernst (5.3) conduce a una solución de las ecuaciones 
de Einstein-Maxwell (<D£, A£[3>¡;],A<) dadas por 

KW«] = A„[O n ], (5.7b) 
JT^]f-(e 1<b - - 1) 

A ' ~ ^.(í-^-a+yS)' (5Jc) 

interpretable como la superposición de un disco con borde interno central alrededor de un agujero 
negro de Schwarzschild cargada. Ahora bien, usando la tétrada ortonormal (2.35) se pueden escribir 
el tensor energía-momentum Q„¡, en la forma canónica como 

Qab = e c V a V b + p c W a W b , (5.8) 

donde e' y p c son, respectivamente la densidad de energía y la presión del disco y el superíndice "c" 
denota que corresponde a la superposición cargada. En términos de estas cantidades se puede escribir 



^ = 4¿W-Kl«® _ r0 cJ ^ (5 9a) 

p , = 4e^- A 5 [v - D r®^0^, (5.9b) 

donde, como arriba, todas las cantidades son evaluadas en z = 0 + . En términos de esta misma tétrada, 
como se sabe del (Capítulo 2), la densidad superficial de corriente en el disco, I a = [F" z ], se puede 
escribir en la forma canónica como 

h = o-V„ + jW a . (5.10) 

Así, explícitamente, la densidad de energía, la presión azimuthal la densidad newtoniana efectiva y la 
densidad de carga eléctica son, respectivamente: 

£ » - 084" + 1)4 £ "' (5 Ua) 

Pat+p? ...... 

p - = -wrw e - (5 - llb) 



(#r+i) 3 



A = rar, M 3 ^ ( 5 ' llc ) 



a"„ = ; e„, (5. lid) 



con p'i = £,', + p¡¡, siendo 



4aE n (x) f o- 2 (2»-2)!g„(x) | 
o" +1 x 2 " +1 6XP \ 2 2 "" 1 a 2 "x 4 " /' 
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Cuando /} = 1 las expresiones anteriores se reducen, como debe ser, a las expresiones (4.79); las 
correspondientes a la densidad de energía, la presión y densidad newtoniana efectiva de la superpo- 
sición en ausencia de carga eléctrica. Ahora, como e > 0 y f > 1, e c y p c serán siempre positivos. 
Consecuentemente, el tensor energía-momentum estará en completo acuerdo con las condiciones de 
energía fuerte y débil. De otra parte, con el fin de satisfacer la condición de energía dominante se 
requiere que p° < e c , lo cual implica que f > (1 + -y'l + 8/J 2 )/(2yS) y que el radio del borde interno 
a > m(4/3 2 + 2 Vi + 8yS 2 + 2) 1/2 /(2/3). 

El tensor energía-momentum puede ser interpretado como la supeiposición de dos fluidos carga- 
dos contrarrotantes. En orden a estudiar esto, se escribe Q" b como [22] 

Q" b = e c + U'l V\ + siU'lU b _, (5.12) 

donde 

c c \p 2 ai+P)(K 2 -i-w) 



<#r+i) 4 



\, (5.13) 



son las densidades energía de los dos fluidos contrarrotantes. Las velocidades de contarrotación son 
dadas por [22] 

V a + UW" 

V% = (5.14) 

VT^TT 2 

donde 

p_ = g+g 
e y3(í 2 -l) 

es la velocidad tangencial de contrarrotación. Así, se tienen dos fluidos de polvo con densidades 
superficiales de energías iguales. Ahora, como £ ± > 0, los dos discos de polvo están en completo 
acuerdo con todas las condiciones de energía. 

Ahora, con el fin de mostrar un análisis del comportamiento de las densidades superficial de 
energía de la superposición, se graficará la densidad superficial de energía adimensional é„ = as„ 
en función de la coordenada radial adimensional f = r/a. Así, en la Figura 5.1, Figura 5.2 y en la 
Figura 5.3 se gráfica s„ en función de f para los tres primeros miembros de la familia de soluciones 
correspondientes a la superposición, con n = 1, 2 y 3 respectivamente, para diferentes valores de los 
parámetros m y a„ = a/a". Así, para cada valor de n, se toma m = 2a/3/[4fi 2 + 2-^1 + 8/3 2 + 2] 1/2 , 
primera figura de arriba hacia abajo, luego m = 2a(¡/2[4(¡ 2 + 2-^1+ 8/3 2 + 2] 1 ' 2 , y finalmente m = 
2a/3/3[4y3 2 + 2 Vi + 8/3 2 + 2] 1/2 . En cada caso se gráfica s{f) para 1 < f < 3 con diferentes valores del 
parámetro a„ = a/a". Primero se toma ái = 1, curvas próximas al borde del disco, luego a 2 = 2. Para 
cada valor de á se toman diferentes valores de /?: 0.1, 0.2, . . . , 0.9 y, finalmente, fi = 1, el caso vacío, 
la curva con menor máximo en cada gráfica. Como se puede ver, el comportamiento de la densidad 
superficial de energía correspondiente a la superposición cargada es similar al comportamiento de 
la densidad de la superposición sin carga. Es decir, en todos los casos la densidad superficial de 
energía es siempre positiva, con un máximo cerca del borde de los discos, y luego decrece rápidamente 
conforme f decrece. Se puede ver que, para un valor fijo de n y m, conforme los valores de á„ crecen, el 
valor del máximo disminuye y se desplaza hacia donde crecen los valores de f. De otra parte, para un 
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valor fijo de á„ y m el máximo crece y se desplaza en dirección de crecimiento de r, conforme n crece. 
Mientras que para oyó, fijos, los máximos crecen como m decrece. Ahora, a medida que la carga 
eléctrica aumenta, /} disminuye, el máximo de la función aumenta, describiendo el comportamiento 
asintótico, solución de vacío, cuando (3 = 1. En la Figura 5.1, Figura 5.2 y en la Figura 5.3 se gráfica 
la densidad de carga eléctrica -<x„, como se puede ver su comportamiento es similar al de la densidad 
de energía, por tanto toda la descripción anterior es también valida para la densidad de carga eléctrica. 

Las soluciones obtenidas representan la superposición cargada de un disco infinito delgado anular 
alrededor de un agujero negro de Schwarzschild. La masa total del disco es finita y el tensor energía- 
momentum está de acuerdo con todas las condiciones de energía. Además los estados asintóticos 
de la superposición cargada se reducen a la supeiposición sin carga eléctrica descrita en el capítulo 
anterior. La solución puede interpretarse como la descripción de un disco constituido por fluidos 
provistos de carga eléctrica contrarrotando alrededor de un agujero negro y que, están también en 
completo acuerdo con las condiciones de energía. 

Las soluciones anteriores presentan algunas propiedades interesantes. Primero, las fuentes de ma- 
terial asociado presenta un comportamiento muy razonable y su tensor energía-momentum está en 
completo acuerdo con las condiciones de energía. Además, su relativa simplicidad cuando son expre- 
sadas en términos de coordenadas esferoidales prolatas y oblatas, permite el fácil estudio de aspectos 
dinámicos, como por ejemplo el movimiento de las partículas cargadas dentro y fuera del disco, la 
estabilidad de las órbitas y la posible existencia de singularidades. 
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Figura 5.1: Densidad superficial de energía e = as en función de f = r/a para el primer modelo, 
n = 1, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = lapUAp 1 + 2 Vi + 8/3 2 + 2] 1/2 , primera figura de arriba hacia abajo, luego 
m = 2a/3/2[4/3 2 + 2^1 + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica s(f) para 1 < f < 3 con diferentes valores del parámetro &„ = a /a". Primero se toma ¿1 = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de á se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.9 y, finalmente, /? = 1, el caso vacío, la curva con menor amplitud en cada gráfica. 
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Figura 5.2: Densidad superficial de energía e = ae en función de f = r/a para el segundo modelo, 
n = 2, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = primera figura de arriba hacia abajo, luego 

m = 2a/3/2[4/3 2 + 2^1 + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica s(f) para 1 < f < 3 con diferentes valores del parámetro a„ = a /a". Primero se toma ¿1 = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de á se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.9 y, finalmente, /? = 1, el caso vacío, la curva con menor amplitud en cada gráfica. 
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Figura 5.3: Densidad superficial de energía e = ae en función de f = r/a para el tercer modelo, 
n = 3, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = 2af3/[4/3 2 + 2 Vi + 8/3 2 + 2] 1/2 , primera figura de arriba hacia abajo, luego 
m = 2ay3/2[4y3 2 + 2 Vi + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica s(f) para 1 < f < 3 con diferentes valores del parámetro &„ = a /a". Primero se toma ¿1 = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de a se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.9 y, finalmente, /? = 1, el caso vacío, la curva con menor máximo en cada gráfica. 
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Figura 5.4: Densidad superficial de carga -a = -acr en función de f = r/a para el primer modelo, 
n = 1, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = 2af3/[4/3 2 + 2 Vi + 8/3 2 + 2] 1/2 , primera figura de arriba hacia abajo, luego 
m = 2a/3/2[4/3 2 + 2^1 + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica <x(F) para 1 < f < 3 con diferentes valores del parámetro á„ = a /a". Primero se toma a¡ = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de a se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.8 y, finalmente, /? = 0.9, la curva con menor máximo en cada gráfica. 
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Figura 5.5: Densidad superficial de carga -cr = -acr en función de f = r/a para el segundo modelo, 
n = 2, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = 2af3/[4/3 2 + 2 Vi + 8/3 2 + 2] 1/2 , primera figura de arriba hacia abajo, luego 
m = 2a^/2[4y3 2 + 2 Vi + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica <r(F) para 1 < f < 3 con diferentes valores del parámetro &„ = a /a". Primero se toma 0-1 = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de a se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.8 y, finalmente, /? = 0.9, la curva con menor máximo en cada gráfica. 
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Figura 5.6: Densidad superficial de carga -& = -acr en función de f = r/a para el tercer modelo, 
n = 3, de la superposición cargada de discos con borde interno central de radio a alrededor de un 
agujero negro de masa m = 2af3/[4/3 2 + 2 Vi + 8/3 2 + 2] 1/2 , primera figura de arriba hacia abajo, luego 
m = 2a/3/2[4/3 2 + 2^1 + 8/3 2 + 2] 1/2 , y finalmente m = 2a/3/3[4/3 2 + 2^1 + 8yS 2 + 2] 1/2 . En cada caso se 
gráfica <r(F) para 1 < f < 3 con diferentes valores del parámetro á„ = a /a". Primero se toma a¡ = 1, 
curvas próximas al borde del disco, luego ¿5-2 = 2. Para cada valor de a se toman diferentes valores de 
j8: 0.1, 0.2, . . . , 0.8 y, finalmente, /? = 0.9, la curva con menor máximo en cada gráfica. 
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Conclusiones Generales 



e presentó un análisis detallado del tensor energía-momentum de un disco delgado estacionario 



axialmente simétrico. Se consideró el modelo más general de discos delgado con esfuerzo radial 
y presión radial diferentes de cero. Así, en este trabajo de investigación se presenta por primera vez 
el análisis completo de todos los posibles casos de tensor energía-momentum que pueden obtenerse a 
partir de soluciones estacionarias axialmente simétricas de las ecuaciones de Einstein-Maxwell en el 
vacío. 

Se expresó el tensor energía-momentum superficial (SEMT) en términos de la tétrada LSO y así 
se obtuvieron expresiones explícitas para las variables dinámicas que caracterizan el SEMT. Esto es, 
expresiones para el vector velocidad de los discos, así como para la densidad superficial de energía, 
esfuerzos principales y función flujo de calor que son válidas para toda clase de discos. De igual modo 
se presentó explícitamente una expresión para la densidad superficial de corriente. Las expresiones 
obtenidas aquí, son entonces la generalización de las correspondiente expresiones para los casos está- 
ticos que fueron presentados en [22] y también generaliza los resultados obtenidos en [29] para discos 
delgados estacionarios. 

También se presentó la generalización estacionaria del Modelo de Contrarrotación (CRM), el cual 
fue desarrollado previamente sólo para el caso estático en [22, 21]. Así entonces, ahora es posible 
obtener expresiones explícitas para todas las cantidades pertinentes al CRM presentes cuando existe 
flujo de calor y para esfuerzos principales arbitrarios. Se consideraron velocidades de contrarrotación 
iguales y opuestas y contrarrotación a lo largo de electrogeodésicas. 

Se obtuvo una ligadura general sobre las velocidades de contrarrotación, necesaria para expresar 
el tensor energía-momentum de los discos como la superposición de dos fluidos perfectos contrarro- 
tantes. La ligadura obtenida es la generalización de la obtenida en [29], para discos sin presión radial y 
sin flujo de calor, donde sólo se consideró fluidos contrarrotantes circulando a lo largo de geodésicas. 

Se encontró también que, en general, no existe posibilidad de tomar las dos velocidades tangencia- 
les de contrarrotación iguales y opuestas ni suponer que los fluidos circulen a lo largo de geodésicas 
(o electrogeodésicas). Entonces, para discos construidos a partir de métricas estacionarias axialmente 
simétricas, las velocidades de contrarrotación no están determinadas de manera única y existe la li- 
bertad que pueda ser usado para obtenerse diferentes modelos de contrarrotación para el mismo disco. 
El desarrollo teórico que condujo a estos resultado constituye el fundamento central del trabajo [27]. 

Se presentó una familia infinita de discos de polvo cargado axialmente simétricos de extensión 
finita con densidades superficiales de energía y carga bien comportadas. Los discos fueron obtenidos 
a partir de soluciones de las ecuaciones de Einstein-Maxwell para un espacio-tiempo conformestático. 
Con el fin de obtener las soluciones, se supuso una dependencia funcional entre el potencial eléctrico 
y el potencial gravitacional y, entre esta y una función auxiliar. Las soluciones fueron expresadas 
entonces en términos de una solución de la ecuación de Laplace correspondiente a la familia de bien 
comportados discos newtonianos delgados de radios finitos, los discos generalizados de Kalnajs [30]. 
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Estos discos delgados tienen una densidad de carga que es igual, salvo por un signo, a su densidad 
de energía, y así constituyen un ejemplo de los comúnmente llamados configuraciones de equilibrio de 
"polvo eléctricamente contrapuesto". La densidad es siempre positiva y bien comportada, anulándose 
en el borde del disco. Además, conforme el valor de m crece, la densidad de energía se concentra 
más en el centro de los discos, presentando un máximo en r = 0 para todos los valores de Sin 
embargo, para los primeros dos modelos con m = 1 y m = 2, para valores pequeños de \k\ la densidad 
de energía presenta un máximo cerca al borde del disco mientras que para valores grandes de el 
máximo se encuentra en el centro del disco. Además, como la densidad de energía de los discos es 
siempre positiva y los discos están hechos de polvo , todos los modelos están en completo acuerdo 
con todas las condiciones de energía. Un hecho de particular relevancia en el estudio de modelos de 
discos relativistas. En efecto, como ya fue mencionado, muchos de los modelos de discos relativistas 
delgados que han sido estudiados en la literatura no están en completo acuerdo con estas condiciones. 

De otra parte, como se puede ver de las ecuaciones en Sec. 3.2, el procedimiento aquí presentado 
puede ser aplicado no sólo a espaciotiempos conformestáticos sino que también puede ser usado para 
obtener soluciones no axialmente simétricas de las ecuaciones de Einstein-Maxwell en el vacío. Por 
consecuencia, estos modelos de discos delgados aquí presentados se pueden generalizar considerando 
para la función auxiliar U soluciones de la ecuación de Laplace sin la imposición previa de simetría 
axial. Así, se está trabajando en esta dirección y los resultados se presentarán en un próximo trabajo. 
Análogamnete, la generalización a espaciotiempos conformestacionarios con campos magnéticos está 
en consideración. 

Se presentó una familia infinita de soluciones de las ecuaciones Einstein en el vacío asitóticamente 
planas y regulares en todas partes. Estas soluciones describen una familia infinita de discos delgados 
con un borde interno central constituidos por materia de polvo, cuya densidad de energía es siempre 
positiva y bien comportada, de tal modo que su tensor energía-momentum está en completo acuerdo 
con todas las condiciones de energía. Además, aunque los discos son de extensión infinita, todos 
tienen masa finita. Ahora, como las funciones métricas de las soluciones se calcularon explícitamente, 
estas constituyen la primera familia de soluciones exactas explícitas completamente integradas para 
tal clase de fuentes de discos delgados. Además, su relativa simplicidad cuando son expresadas en 
coordenadas esferoidales oblatas, permite estudiar muy fácilmente sus diferentes aspectos dinámicos, 
como el movimiento de partículas dentro y fuera de los discos y la estabilidad de las órbitas. 

Además de su importancia como una nueva familia de soluciones exactas y soluciones explícitas 
de las ecuaciones de Einstein en el vacío, la principal importancia de esta familia de soluciones es que 
puede ser fácilmente superpuesta con la solución de Schwarzschild con el fin de obtener la descripción 
de un sistema binario compuesto por un disco delgado alrededor de un agujero negro central. En 
efecto, la superposición del primer miembro de esta familia con un agujero negro de Swarzschild fue 
obtenida y presentada en [31]. 

Ahora, como en la teoría newtoniana, el potencial gravitacional está dado por la solución del 
problema de valores en la frontera para la ecuación de Laplace, se puede considerar <£„(/■, z; a) como 
una familia de potenciales gravitacionales newtonianos de fuentes discoidales con borde interno, cuya 
densidad de masa newtoniana están dada por 

2aE„(x) 

claramente divergente en el borde de los discos. Por consecuencia, se puede concluir que no existen 
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soluciones en la teoría newtoniana que describan el campo gravitacional de un disco delgado con 
borde interno, mientras que esta clase de fuentes puede ser propiamente descrita por medio soluciones 
regulares y asintóticamente planas, en el marco de la teoría de la gravitación ofrecida por la relatividad 
general. Los fundamentos básicos necesarios para la obtención de esta clase de soluciones fueron 
presentadas en el trabajo [25]. 

Se presentó también un formalismo para superponer discos de polvo con un borde interno central 
alrededor de un agujero negro de Schwarzschild. La superposición se obtuvo usando la solución de 
las ecuaciones de campo de Einstein correspondientes a discos de polvo con un borde interno central 
del trabajo [25]. El formalismo para la superposición fue inspirado en el trabajo de Letelier and 
Oliveira [56]. Las densidades de energía del sistema superpuesto son siempre funciones positivas de la 
coordenada radial, igual a cero en el borde interno central del disco, con máximo en un valor finito de 
los radios y cayendo a cero en el infinito. Los discos tienen masa finita y su tensor energía-momentum 
está de acuerdo con todas las condiciones de energía. La presión azimutal y la densidad newtoniana 
efectiva de los discos tienen comportamiento similar al de la densidad de energía, así que el análisis 
previo es aplicable también a estas cantidades. De otra parte, la velocidad de contrarrotación de las 
partículas es siempre positiva y menor que la velocidad de la luz, cayendo a cero en el infinito. La 
masa del disco cuando el agujero negro está presente es menor que la masa total del disco solo. 

Finalmente, se presentó, lo que constituye hasta donde se sabe, la primera familia de soluciones 
exactas completamente integradas explícitamente de las ecuaciones de Einstein-Maxwell correspon- 
dientes a la superposición de un disco de polvo contrarrotante con un agujero negro. Las soluciones 
obtenidas representan la superposición cargada de un disco infinito delgado anular (un disco con un 
borde interno) alrededor de un agujero negro de Schwarzschild. La masa total del disco es finita y el 
tensor energía-momentum está de acuerdo con todas las condiciones de energía. Además los estados 
asimptóticos de la superposición cargada se reducen a la superposicón de sin carga eléctrica. La so- 
lución puede interpretarse como la descripción de un disco constituido por fluidos provistos de carga 
eléctrica contrarrotando alrededor de un agujero negro y que, están también en completo acuerdo con 
las condiciones de energía. 

Como se puede ver, las soluciones anteriores presentan algunas propiedades interesante. Primero, 
las fuentes de material asociado presenta un comportamiento muy razonable y su tensor energía- 
momentum está en completo acuerdo con las condiciones de energía. Además, su relativa simplicidad 
cuando son expresadas en términos de coordenadas esferoidales prolatas y oblatas, permite el fácil 
estudio de aspectos dinámicos, como por ejemplo el movimiento de las partículas dentro y fuera del 
disco, la estabilidad de las órbitas y la posible existencia de singularidades. Ahora bien, aunque el 
formalismo de los potenciales complejos de Ernst fue usado aquí para obtenerse la superposición 
cargada, también puede ser usada para generar una superposición magnetizada y además una cargada 
y magnetizada. Resultados sobre este particular se presentarán en un próximo trabajo. 

Los principales resultados obtenidos en el marco del desarrollo del presente trabajo de investiga- 
ción, así como trabajos en curso, nacidos desta misma investigación, se listan a continuación: 

1. An Infinite Family of Magnetized Morgan-Morgan Relativistic Thin Disks [34], 

2. An infinite family of relativistic magnetized finite thin disks [35], 

3. Stationary Axially Symmetric Relativistic Thin Disks [26], 
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4. Charged and Electromagnetized Morgan-Morgan-Type Disks [36] , 

5. Modelos Relativistas de Discos Delgados Finitos No Axialmente Simétricos con Campo Eléctri- 
co, Presentación Oral en Segunda Reunión Colombo-Venezolana de Relatividad y Gravitación, 
Armenia, Colombia 28 al 31 de Octubre de 2007, 

6. Discos de Morgan-Morgan Cargados y Electromagnetizados, Presentado en Calidad de Poster 
en el Tercer Taller de Gravitación y Astrofísica Relativista, Homenaje a Luis Herrera-60 Aniver- 
sario, Isla de Coche, Venezuela 8 al 11 de Noviembre de 2006, 

7. Finite Axisymmetric Charged Dust Disks Sourcesfor Conformastatic Spacetimes [23], 

8. Relativistic Static Thin Dust Disks with an Inner Edge: An infinite Family of New Exact Solu- 
tions [25], 

9. Finite Axisymmetric Charged Dust Disks in Conformastatic Spacetimes, Physical Review D 
[27], 

10. Exact Static Axially Symmetric Thin Annular Dust Disks [24 ], 

1 1 . Electrovacuum Static Axially Symmetric Thin Annular Dust Disks, Presentado en Calidad de 
Poster en // Congreso latinoamericano de Física, 20 al 24 de Octubre, Zacatecas, México, 

12. Counterrotating Dust Disk Around a Schwarzschild Black Hole: New Fully Integrated Explicit 
Exact Solution [31], 

13. A Family of Relativistic Charged Thin Disks with an Inner Edge [66], 

14. Counterrotating Charged Dust Disk Around a Schwarzschild Black Hole: New Fully Integrated 
Explicit Exact Solution [65], 

Los tópicos aquí estudiados también sirvieron como punto de partida para los trabajos de grado: 

1. Trabajo de Grado Para Optar el Título de Maestría en Física Construcción de Modelos Re- 
lativistas de Discos Delgados Mediante Soluciones Conformestáticas de las Ecuaciones de 
Einstein-Maxwell, cuya calificación obtenida fue meritoria. 

2. Trabajo de Grado Para Optar el Título de Físico Modelos Relativistas de Discos de Polvo axial- 
mente Simétricos con Borde Interno, 

codirigidas por los autores del presente trabajo de investigación. 
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Apéndice A 

La Ecuación Dinámica y la Ecuación de la 

Electrogeodésica 



A.l. La Ecuación Dinámica 

En este apéndice se obtiene una "ecuación dinámica" asociada a los campos gravitacional y elec- 
tromagnético debido a fuentes concentradas en la hipersuperficie discoidal. El punto de partida es la 
ley de conservación 

r\=0, (A.l) 



jai, _ Q^S^ + iTfj 0 , (A.2a) 

T t = Tfj z ,®(<p) + Clrd " (A-2b) 

donde &((/>) denota la función de Heaviside con &(</> = 0) = =, 

Sustituyendo el tensor energía-momentum (A. 2) en la ecuación de conservación (A.l), se tiene 
que la derivada covariante de T,,/, es 

r% = (G*W)i* + {t£ l-®(0) + t;X^ - ■ (A.3) 

Ahora bien 

T a %+ . b = T'\- ,, = 0, 

EM lz '° EM lZ ' ü ' 

de donde se deduce que 

r* 4 = (Q" b ó(<f>)), b + [Tf u ] <pj,d($). (A.4) 
El primer término del lado derecho de esta ecuación es 

wfsyifh = Q ,b ., b m + <r b s'(m.b, (a.5) 

Entoces, las ecuaciones (A.l) con 

T ab , b = { &'";„ + [T£\ <P¿} m + Q'" '<t>j,8\</>) = 0, (A.6) 

Así se tiene que 

= 0, (A.7a) 

Q"",,, + [2£]*» = 0. (A.7b) 
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Por lo que se puede obtener 

+ [T£\ = 0. (A.8) 

donde se han usado la definición de los vectores normales a la superficie ¡p¡, = 5 z b . De otra parte, se 
sabe que 

T"¿ = F" C F ZC - l -g" z F cd F«< (A.9) 

entonces se tiene que el salto en el tensor energía-momentum correspondiente al tensor electromag- 
nético a través de la superficie del disco es 

[T a ¿] = í F "cl F '~ C + F "c H - \<t [ F a [F'"] + [F ci \ F"'} . (A. 10) 

o, equivalentemente, 

[T EU 3 = [F ac ] F zc + F ac [F zí ] - I¿¡ { [F c ¿ F c " + F cá [f"'\ } . (A. 1 1 ) 

Ahora, debido a la propiedad F c ¿ [f c<í J = [F«¡] F cd se tiene 

IT» 3 = [KA F* + F ac [F a ] - l -ólF"' [F cd ] . (A. 1 2) 

Es muy fácil darse cuenta que F zc [F ac ] = 6 z a F zc [A c ,] y que [F cí ¡] F cd = 2F ZC [A cz ]. Así, se obtiene 
finalmente que para el salto del tensor electromagnético 

[T'« a = F ac [F zc ] . (A. 13) 

Por tanto, se tiene de (A. 8) que 

Q a "., b = F ac [F cz ]. (A. 14) 

La ecuación (A. 14) puede ser renombrada como una "ecuación dinámica" asociada al campo gra- 
vitacional y gravitacional debido a fuentes concentradas en la hipersuperficie discoidal. El punto de 
partida para su obtención es la ley de conservación T ., = 0. 

A.2. La Ecuación de la Electrogeodésica 

Ahora se supone que el disco está constituido por la superposición de dos fluidos contrarrotantes, 
lo cual es, como ya se sabe, equivalente a escribir las relaciones 

S a b = S.l + S„ b _, (A. 15a) 

j" = i t + i-, (A. 15b) 

donde 

Sf = (c ± + P± )UIUÍ + P± h" b , (A.16) 
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Con, S" b = e K ~®Q lh , de la ecuación (A. 16) se tiene que 

Qf = (€ ± + p ± )u:u b ± +p ± h"\ (A.17) 

donde k = e®~ K K. Después de un cálculo medio tedioso y algo aburrido, podemos obtener la derivada 
covariante de Q¿, la cual es: 

Qj. b = (g ± + p ± )U ± a ., b U b ± + p ± ,„ h¿> + p ± h a ». b , (A.18) 

así que, después de calcular la forma explícita de h Jr. b se tiene 

Qj\t = + P ± )V b ± V L ± g bcf¡ + P±JI + p ± 6\ (A - O),,.. (A.19) 

Ahora, de la ecuación (A. 14) se tiene que 

;» = ■"*£. ( A - 2 °) 

así, por reemplazar (A.19) en (A.20) se tiene 

- l -(e ± + p ± )U»Ulg bcjl + p ± , a + p ± <5 r ( ,(A - O),,. = F a ,cr ± U¿ (A.21) 

donde, como antes, k = e®~ A K. De (A.21) esperaríamos cuatro ecuaciones, una para cada valor de a. 
Sin embargo sólo se tiene una, dado que, excepto para x 2 = r, se tienen ecuciones triviales. Así se 
tiene 

- ife + p ± )U':u<g bc , + p ±J = F rc o- ± Ul (A.22) 

Ahora, A,. = 0, de modo que F n = A, ,-. Consecuentemente, con la sustitución de esta cantidad en la 
ecuación anterior se obtiene 

- -(e ± + p ± )U ± b U^g bL , r + p ± , = A„.o- ± u:. (A.23) 
Ahora se desarrolla cada una de las sumatorias en (A.23). La primera: 

U b u:g bíV . = [/ ± °[/ ± %o,,- + 2U°U¡g 0U + U¿ U>g n , r , (A.24) 
El término de la derecha en (A.23) es, explícitamente, 

A c , r U¿ = A 0ir U° + A hr U¡, (A.25) 
luego, por sustituir (A.24) y (A.25) en (A.23) se obtiene finalmente 

\(U*) 2 ( gaa , r + 2g m , r O-± + gn,M\){e ± + p±) - P ± , = -o- ± U° ± (A,, r + A 9 , r Q ± ) , (A.26) 

donde se ha llamado Q ± = U l ± IU a ± , y como se sabe A„ = (A„A 9 ,0,Q). La ecuación (A.26) es la 
ecuación de la electrogeodésica para un disco delgado estacionario constituido por la superposición de 
dos fluidos contrarrotantes cargados eléctricamente, con presión radial. Costituye una generalización 
de la ecuaciones en la electrogeodésica presentada en [21]. 
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